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Abstract 

This  thesis  examined  the  effects  of  system  configura¬ 
tion  and  transmitted  waveform  as  applied  to  an  airborne 
bistatic  radar.  The  study  concentrates  on  the  limiting 
effects  of  the  transmitted  signal  parameters,  the  signal  to 
clutter  ratio,  and  the  doppler  spread  of  the  receiver  beam. 

The  analysis  was  performed  by  constructing  a  computer 
model  to  graphically  display  the  results  of  the  ai^licable 
equations.  The  equations  developed  in  this  study  arei  The 
power  limited  region  equation,  the  isorange  equation,  the 
isodoppler  equation,  and  the  signal  to  clutter  ratio.  These 
equations  and  those  of  ellipsoidal  geometry  were  used  to 
develop  the  bistatic  radar  model. 


AIRBORNE  BISTATIC  RADAR 
LIMITATIONS  Al>rD  SAMPLE  CALCULATIONS 

I .  Introduction  to  Bistatic  Radar 
Airborne  Bistatic  Radar 

An  Airborne  Bistatic  Radar  System  is  an  electro-magneti 
detection  device,  composed  of  separate  transmitter  and 
receiver  platforms,  used  to  find  range  and  velocity  charac~ 
teristics  of  objects  by  measuring  time  of  flight  and  doppler 
shift  of  the  transmitted  signal.  The  receiver  is  alvays 
mounted  on  an  airborne  platform. 

Assumptions 

A  basic  assumption  of  this  thesis  requires  that  the 
transmitter  and  receiver  are  syntd^ionized  in  some  fashion. 
Further,  it  is  assumed  that  the  receiver  processor  has  the 
ability  to  detect  the  signal  and  do  whatever  processing  is 
required.  These  assumptions  were  set  by  the  Tactical  Sur¬ 
veillance  Division,  Rome  Air  Development  Center,  Griffiss 
AFB,  NY.  This  thesis  concentrates  on  the  limiting  effects 
of  the  transmitted  signal  parameters,  the  signal  to  clutter 
ratio,  and  the  doppler  spread  of  the  receiver  beam.  The 
signal  detection  and  estimation  problem  as  veil  as  tne  syn¬ 
chronization  problem  are  under  investigation  by  others. 
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The  Problem 


The  airborne  receiver  platform  motion  is  not  constrained 
to  any  velocity  or  motion  with  respect  to  a  transmitter. 
Therefore,  the  received  signal  is  subject  to  the  inherent 
characteristics  resulting  from  bodies  with  independent 
motion j  in  particular,  the  doppler  shift  of  any  signal  whicdi 
would  travel  between  the  transmitter  and  the  receiver,  or 
between  the  transmitter,  a  target  and  the  receiver  platform. 
This  thesis  outlines  the  characteristics  of  a  signal  which 
limit  the  surveillance  area  of  the  receiver  platform.  A  com¬ 
puter  model  has  been  developed  to  incorporate  the  concepts 
outlined  in  this  study.  The  model  (see  Appendix  A)  provides 
an  analysis  tod  for  bistatic  radar  development. 

Developments  in  the  radar  field  have  moved  rapidly  for¬ 
ward  with  the  marriage  of  the  digital  computer  to  the  radar 
detection  arena.  The  terminology  in  a  bistatic  radar  system 
is  built  upon  a  traditional  base.  However,  new  terminology 
must  be  introduced  to  explain  the  change  in  the  problem  and 
system  configuration.  The  terminology,  limitations  and  major 
components  of  the  bistacic  system  are  outlined  in  chapter  II. 
The  configuration  of  the  platforms  and  characteristics  of 
the  transmitted  waveform  are  discussed  as  well  as  some  advan¬ 
tages  of  a  bistatic  system.  A  list  of  current  research  is 
also  provided. 

ji^nalysis  of  the  Problem 

~  The  defining  param^ers~^.3  'equations  “are '^e^l oped 
chapter  III.  In  particular,  tlie  range,  geometry  and  system 
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performance  equations  are  developed  to  demonstrate  the  sys¬ 
tem  configuration  as  well  as  the  ability  of  the  system  to 
perform. 

Limiting  Regions .  The  limiting  regions  of  the  bistatic 
radar  are  defined  by  two  equations.  The  power  and  the  iso¬ 
range  equations  are  developed  to  explain  the  viewing  area  of 
a  bistatic  radar  and  its  unique  blind  zone.  For  example, 
there  is  a  minimum  and  a  maximum  range  associated  with  the 
isorange  equation  which  is  independent  of  the  power  equation. 

Doppler  Effects.  The  bistatic  doppler  equation  is  not 
the  same  as  the  monostatic  radar  doppler  equation.  The 
motio’^  of  the  three  objects  (receiver,  transmitter,  and  tar¬ 
get;  must  be  considered  in  their  relative  radial  velocities 
with  respect  to  each  other.  The  doppler  equation  must  also 
account  for  the  secondary  doppler  shift  (and  relative  change 
in  wavelength)  as  the  signal  travels  from  the  transmitter, 
strikes  the  target,  and  is  reradiated  to  the  receiver  plat¬ 
form. 

Signal  to  Clutter  Ratio.  The  signal  to  clutter  power 
ratio  (SCR)  is  defined  and  parameters  of  interest  are 
explained.  The  clutter  cell  area  is  defined  and  its  affect 
on  the  SCR  is  shown.  A  surface  cross  section  per  unit  area 
is  modeled  from  previous  work  and  incorporated  into  the  model 

Finally  the  doppler  spread  of  the  clutter  cell  is 
developed.  The  equations  for  the  doppler  spread  show  that 
the..bis.tatic.  sy^stem  is  dynamic  and  dependent  upon  the  motion 
and  beam  shapes  of  the  transmitter  receiver  pair. 
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Demonstration  of  the  Model 

The  results  of  the  analysis  were  put  into  an  inter¬ 
active  computer  model  which  was  exercised  to  demonstrate 
some  Of  the  features  outlined  in  chapters  II  and  III. 

Chapter  IV  contains  the  results  of  the  simulation.  The 
waveform  limited  regions  are  shown  to  be  dependent  upon  the 
pulse  repetition  frequency  (prf),  the  pulse  width  and  the  prf 
multiple.  The  effects  of  these  variables  are  shown  in 
graphical  form.  Doppler  contour  maps  of  the  system  are 
shown  to  explain  the  effects  of  the  transmitter  receiver 
target  motion.  The  maps  demonstrate  the  effects  of.  the 
relative  velocities  of  the  transmitter  and  receiver  on  the 
ground  plane. 

Signal  to  clutter  tables  and  graphs  are  provided..  The 
effect  of  changes  in  target  radar  cross  section  on  the  SCR 
is  demonstrated.  The  doppler  spread  of  the  receiver  beam  is 
also  shown  in  graphical  form.  It  is  demonstrated  by  the 
graphs  that  the  effects  of  the  platform  motion  as  well  as 
the  receiver  look  angles  have  predictable  results.  The 
doppler  spread  is  shown  to  be  a  well  behaved  phen(»nenon. 
Finally,  a  sample  air  to  air  senario  is  provided  to  show  the 
dynamic  nature  of  the  system. 

Findings  and  Conclusions 

A  summary  of  the  detect ioi  regions  and  the  minimizaticm 
or  maximization  or  the  detection  regions  is  provided  in 
chapter  V»  Several  reconwenda t ions- for- foi low -cn -analysis — 
are  provided. 
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II.  Bistatic  Radar 


Background 

The  idea  of  a  passive  radar  excites  the  imagination  of 
the  Air  Force  information  gathering  community.  If  it  vere 
possible  to  build  a  radar  system  which  used  separate  trans¬ 
mitter  and  receiver  locations,  extra  protection  would  be 
afforded  to  the  radar  receiver  because  the  radar  receiver 
would  not  broadcast  its  presence.  Further,  the  passive 
radar  receiver  would  be  less  vulnerable  to  classic  anti-radar 
teclmir  ’.es  such  as  anti -radar  missiles  or  electronic  counter 
measures.  A  passive  radar  system  would  require  the  presence 
of  uoable  host  signal.  The  host  signal  would  be  coopera¬ 
tive  if  the  transmittex-  and  receiver  were  working  together, 

A  non-cooperative  host  would  be  a  transmitter  whose  broad¬ 
cast  energy  were  being  used  without  the  transmitter’s  know¬ 
ledge. 

Several  studies  have  oeen  commissioned  to  explore 
methods  of  realising  a  passive  bistatic  radar  system.  For 
example,  the  study  ox  bistatic  radar  cross  sections  (5j  6j 
8),  separate  host  recei\<-r  synchronization  methods,  1|  3}  4j 
10),  and  a  synthetic  ape.  ture  mapping  demonstration  (1)  are 
available.  XonTech  explored  the  possibility  of  using  a 
non-cooperative  host  acquisition  technique  and  delineated 
the  receiver  criteria  xor  using  such  a  separate  host 
receiver  system  (4j21). 


Historic  BacKqround  of  Radar  Systems .  The  historic 
background  of  radar  systems  is  flecked  vith  use  of  iriulti- 
static  systems .  The  early  systems  developed  in  Germany  and 
England  were  bistatic.  An  array  of  receive  antennas  were 
placed  at  known  distances  from  the  transmitter  antenna.  The 
radar  returns  were  monitored  and  useful  information  was 
obtained  by  triangulation  between  several  receivers.  The 
use  of  fence  radar  systems  employ  this  technique.  Some  of 
the  Early  Warning  nets  used  this  technique  (12:36-2), 

During  the  Second  World  War,  technological  breakthroughs 
eliminated  the  need  for  separate  receiver  and  transmitter 
antenna.  As  knowledge  of  radar  waves  increased,  techniques 
were  developed  for  protecting  the  receiver  from  the  trans¬ 
mitter  antenna  power  during  transmission.  Techniques  used 
included  the  use  of  recirculators  and  duplexers.  The  elimi¬ 
nation  of  the  need  for  the  separate  transmitter  and  receiver 
antennas  led  to  the  development  of  the  monostatic  radar. 

These  systems  were  attractive  due  to  the  improvement  in  cal¬ 
culation  of  range  and  doppler  information.  The  classical 
radar  range  and  doppler  equations  (lit 4)  provided  useful 
range  and  velocity  values.  The  radars  were  useful  in  devel¬ 
oping  air  control  radar  systems.  As  more  radar  propagation 
knowledge  was  gained,  it  became  feasible  to  mount  a  radar  on 
an  airborne  platform.  The  radar  was  used  to  provide  naviga¬ 
tion  and  surveillance  information  (14j6). 

Historically,  the  processing  of  the  radar  returns  was 
done  with  analog  computing  devices.  The  marriage  of  the 
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dioitvii  processor  to  tlio  radar  family  has  greatly  expanded 
radar  systesn  capabilities,  Th/s  growing  level  of  sophistica¬ 
tion  of  the  Electronic  Counter  Measure  (ECM)  capabilities 
has  generated  a  new  interest  in  the  Bistatic  Radar  System. 

The  Bistatic  Radar  returns  are  subject  to  the  same  type  of 
counter  measures,  but  the  ECM,  for  effective  bistatic  black¬ 
out  ,  would  require  a  system  much  broader  in  power  and  techni¬ 
cal  complexity  (9:24) . 

Probl em 

An  important  aspect  of  the  bistatic  radar  problem  is 
the  interpretation  of  the  data  which  the  system  delivers. 

The  bistatic  system  does  not  use  the  classical  radar  range 
equation  (121557).  Additionally,  the  system  will  perceive 
different  doppler  shift  (change  in  transmitted  signal  fre¬ 
quency)  effects  due  to  the  movement  between  the  transmitter, 
the  target,  and  the  receiver  (15*25).  The  classical  model 
has  a  doppler  shift  which  is  twice  the  velocity  divided  by 
the  transmitted  wave  length  (2V/L)  (14*592).  The  doppler 
shift  in  bistatic  model  is  due  to  the  combined  effects  of 
the  radial  velocity  between  the  transmitter  and  target 
divided  by  the  transmitted  wavelength  (VoA»o)  and  the  radial 
velocity  between  the  target  emd  receiver  divided  by  the  wave¬ 
length  of  the  reflected  frequency  <V1 AN) .  The  system -is 
best  understood  using  the  concepts  of  isoranges  and  isodops. 
The  isorange  curve  gives  regions  of  constant  range  and  is  the 
solution  of  the  bistatic  radar  range  equation  for  a  specific 


^ et  of  input  values.  The  isodop  is  a  curve  of  constant 
velocity  which  results  from  the  bistatic  doppler  equation. 

The  purpose  of  this  study  is  to  build  an  analytic  com¬ 
puter  program  model  to  study  the  airborne  bistatic  radar 
case.  The  model  win  be  used  to  develop  isorange  and  isodop 
maps  to  explore  the  possible  airborne  bistatic  system  con¬ 
figurations  (e.g.,  Air  to  Air,  Air  to  Ground). 

The  bistatic  radar  problem  can  be  solved  for  many  con¬ 
figurations  for  the  transmitter  receiver  pair.  The  model 
developed  in  this  study  is  general  in  scope  but  is  limited 
to  pulse  doppler  radar  waveforms.  The  general  case  is  devel¬ 
oped  but  the  specific  case  study  is  limited  to  a  High  PRF 
(25,000  Hz),  Multiple  PRF  (4),  S-Band  airborne  receiver. 

The  receiver  is  assumed  to  be  capable  of  acquiring  and  proc¬ 
essing  the  host  transmitter  waveform.  The  actual  method  for 
acquiring  the  signal  is  not  a  part  of  this  study.  This 
study  concentrates  on  limitations  imposed  by  the  signal 
characteristics . 

Definition  of  Key  Terminology 

The  airborne  bistatic  radar  system  is  composed  of  sev¬ 
eral  components  of  the  classic  monostatic  radar  system. 

The  configuration  of  these  components,  and  more  specifically, 
the  roll  of  these  components  is  the  important  issue  in  under¬ 
standing  the  airborne  bistatic  case.  The  signal  (waveform) 
properties  of  the  host  signal  transmitting  system  determine 
the“  "pe rTOr mance"  OT'the  “radaiT'sy s tew;  The  Teeei vef  ’  s “acqu i - 
sition  of  the  signal  requires  techniques  which  vary 
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depending  upon  "Whether  or  not  the  system  is  cooperative  or 
non-cooperative.  For  example,  a  cooperative  system  would  be 
in  constant  contact  (e.g,,  a  secure  data  link)  to  assure 
beam  synchronization  while  a  non-cooperati"ve  system  would 
most  likely  be  constantly  chasing  the  source  beam  in  the 
"expected*'  trace  region.  The  bi.  tatic  system  is  sensitive 
to  the  power  limitations  imposed  by  the  classic  radar  range 
equation.  Assuming  the  receiver  acquires  the  signal,  the 
signal  characteristics  have  a  direct  effect  on  performance 
(range) . 

Ma 1or  Components .  The  major  components  of  the  Bistatic 
system  are  the  transmitter,  target,  and  receiver.  The 
bistatic  receiver  is  separate  from  the  transmitter  by  defini¬ 
tion  (15tl).  The  role  of  each  component  is  important  in 
quantifying  system  performance. 

The  airborne  receiver  is  the  key  element  of  the  bistatic 
system.  It  is  responsible  for  acquiring  the  signals  emitted 
by  the  transmitter .  The  emitted  signal  is  acquired  on  two 
paths j  A  direct  path  between  the  transmitter  and  receiver 
and  the  bistatic  reflected  path  from  the  transmitter  to  the 
target  to  the  receiver.  The  airborne  receiver  platform  is 
assumed  to  be  a  surveillance  -  listening  -  system  which  flys 
a  constant  elevation  orbit  in  a  sanctuary  zone  a  considerable 
distance  (100-200  Kilometers)  from  the  transmitter. 

The  transmitter  is  any  radar  (airborne  or  ground  based) 
which  emits  a  waveform  that  can  be  processed  by  the  receiver. 
The  transmitter  may  be  stationary  or  mobile. 
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The  target  is  any  object  that  reflects  the  transmitted 
energy  and  produces  a  bistatic  reflection  vhich  is  monitored 
by  the  receiver.  The  target  has  characteristics  of  cross 
section  and  velocity  which  affect  the  bistatic  radar  equa¬ 
tions  as  a  function  of  the  radial  distance  and  radial  veloc¬ 
ities  . 

Signal  Properties .  The  bistatic  system  is  capable  of 
operating  using  cooperative  and  non-cooperative  host  trans¬ 
mitters.  Accordingly,  the  acquisition  and  processing  of  the 
signals  is  a  prime  consideration.  A  PRIORI  Knowledge  of  the 
signals  increases  the  accuracy  of  the  information  gleaned 
from  the  output  of  the  receiver.  This  requires  some  practi¬ 
cal  considerations  for  system  performance  characteristics 
which  are  limited  by  bandwidth  and  host  transmitter  power. 

There  are  several  methods  for  signal  acquisition 
posed  in  the  literature.  A  ncm-cooperative  system  requires 
the  receiver  to  search  the  frequency  spectrum  for  a  waveform 
with  the  desired  characteristics  (4i2).  The  waveform  roust  be 
capable  of  being  processed  by  the  receiver  hardware.  For 
example,  a  pulse  doppler  processing  system  may  not  be  able 
to  process  a  ccmtinuous  Wave  Frequency  Modulated  <CW-K1) 
signal.  A  cooperative  system  could  use  flight  path  synchro- 
niaation  (9)  or  secure  data  links  to  periodically  update 
the  system  synchronization  (1). 

Synchronization  of  the  beam  between  the  transmitter  and 
Jthe  receiver  is  important  J 14* The_tjraQS|nitJ:er_^am_ 
systematically  sweeps  Out  a  volume  of  space.  The  bistatic 


10 


signal  reflections  (returns)  available  during  the  sveep  are 
limited  by  the  beam  shape  and  power  density.  A  cooperative 
host-receiver  system  has  the  option  of  a  direct  data  linK  to 
enhance  signal  synchronization  and  information.  Knowledge 
of  the  beam  sweep  pattern  will  allow  the  receiver  to  search 
in  the  expected  return  area.  A  non-cooperative  host  system 
requires  the  same  synchronization.  Because  the  non-coopera¬ 
tive  host  receiver  system  is  not  capable  of  perfect  phase 
synchronization,  the  phase  information  is  not  available  to 
the  receiver  processor.  The  loss  of  phase  limits  the  overall 
system  performance  but  the  loss  of  phase  information  does  not 
preclude  the  use  of  non-cooperative  hosts  in  pulse  doppler 
radar  applications. 

A  method  to  avoid  the  high  cost  of  beam  synchronization 
is  the  use  of  a  multiple  beam  search  antenna  system  (Figure  1) 
An  array  of  beams  searches  the  volume  of  space  where  the 
transmitter  is  providing  illumination  (9t39}  14i577).  Beam 
synchronization  is  limited  to  systems  where  the  ejqpected 
transmitter  sveep  pattern  can  be  replicated.  For  example,  a 
scanning  search  radar  system  has  physical  limitations  in 
antenna  movement  and  pointing  capabilities.  A  leased  array 
which  could  change  wavefront  direction  at  whim  without  con*^ 
sideration  of  antexina  orientation  would  be  a  poor  host. 

Acquiring  the  bistatic  signal  return  does  not  guarantee 
any  meaningful  results  from  the  receiver  data  processing 
system.  17ie  bistatic  return  is^subject  to  mi^tipath  and 
clutter  constraints  as  veil  as  time  of  flight  versus  power 
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Figure  1  Search  With  Multiple  Receiver  Antenna  Beams 

considerations.  Assumptions  about  the  non-cooperative  trans¬ 
mitter  vaveform  limit  the  potential  host  exploitation. 
Further,  a  non-cooperative  host  which  could  reconfigure  its 
waveform  parameters  is  a  poor  host  candidate  for  the  bistatic 
receiver.  The  ability  to  capture  meaningful  data  is  the 
prime  concern  for  the  re<^ver. 

The  extraction  of  information  contained  in  the  bistatic 
signal  return  is  limited  by  the  transmitter,  target  and 
receiver.  The  bistatic  signal  return  contains  information 
in  the  time  domain  which  is  decoded  into  range  information 
of  the  transmitter-target-receiver  family.  Additionally, 


the  bistatic  signal  return  contains  frequency  domain  infor¬ 
mation  which  decodes  into  the  relative  radial  velocities  of 
the  system  components.  The  frequency  domain  information  is  . 
limited  by  the  PRF  and  pulse  width. 

The  transmitter  power  limits  the  detection  region  of  the 
bistatic  system  because  the  receiver  is  assumed  to  be  com¬ 
pletely  passive.  Reliance  cn  a  host  limits  the  range  of  the 
system  under  certain  conditions .  For  example,  Westinghouse 
asserts  that  under  certain  conditions  (Bistatic  baseline 
separation  is  less  than  1.5  times  the  nominal  Monostatic 
target  detection  range)  the  loss  of  coverage  area  relative 
to  the  Monostatic  radar  is  insignificant  (2t3-124).  The 
pas‘=!ive  nature  of  the  receiver  limits  the  radar  range  equa¬ 
tion  variables  controllable  by  the  receiver  (e.g.,  SNR, 
front  end  temperature  K,  antenna  gain,  etc.).  These  vari¬ 
ables  are  usually  hardware  dependent  and  not  dynamically 
adjustable. 

The  Ovals  of  Cassini  are  used  to  describe  the  power 
limitations  of  a  Bistatic  System  (12i3C-41).  The  Ovals  of 
Cassini  describe  the  region  where  the  receiver  is  capable  of 
detecting  useful  information.  Refer  to  Figure  5  and  accom¬ 
panying  text  for  a  complete  development  of  the  Ovals  of 
Cassini. 

The  transmitter  beam  characteristics  affect  the  system. 
The  main  beam  lobe  carries  most  of  the  power  and  is  most 
capable  of  detection.  The  limitations  imposed  by  jamming 
of  the  main  lobe  will  not  be  as  critical  due  to  the  bistatic 
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return  geometry  because  the  bistat.i.c  return  pattern  is  not 
the  same  as  the  monostatic  return  (9i23).  This  study  will 
not  consider  side  lobe  returns. 

Physical  Properties  of  the  Transmitted  Waveform.  Phys 
ical  properties  of  the  transmitted  waveform  are  critical  in 
the  analysis  of  the  bistatic  problem.  A  class  of  pulse  dop 
pier  radar  waveforms  is  examined  in  this  study.  The  class 
is  limited  to  a  generic  waveform  which  is  realizable.  The 
pulse  repetition  frequency,  the  puli?e  width,  and  receiver 
bandwidth  are  critical  design  parameters.  Westinghoase 
asserts t 

These  factors,  as  well  as  legal  frequency 
allocation  requirements,  drive  radar  designs  to 
operating  in  one  restricted  frequency  band  with  an 
instantaneous  bandwidth  only  vide  enough  to  provide 
the  range  cell  length  needed  by  the  systan.  For 
typical  long  range  airspace  surveillance  radars, 
effective  pulse  widths  in  the  vicinity  of  0.5  to  2 
microseconds  limits  the  instantaneous  bandwidth  to 
about  0.5  to  2  MHZ  [2t?-127]]. 

The  Pulse  Repetition  Frequency  (PRF)  is  the  frequency 
at  which  the  transmitter  emits  pulses.  The  VR?  is  the 
reciprocal  of  the  time  between  pulses  (lA).  The  use  of 
multiple  PRFs  (lit  114)  increases  the  range  effectiveness  of 
the  radar. 

Clutter  is  any  radar  return  (monostatic,  bistatic,  or 
multistatic}  which  is  not  desired.  A  large  signal  return 
due  to  background  could  hide  a  possible  target.  Clutter 
echos  greatly  distort  the  information  being  processed  by 
4:be  sysbcm.  One  of  the  secondary.  Concerns  „oX  jthe,bisJ;atic 
radar  range  and  doppler  problem  is  the  need  to  quantify  the 
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clutter  phenomena  of  the  bistatic  system.  A  clutter  thresh¬ 
old  is  usually  established  to  limit  the  unwanted  signal 
returns  on  a  system.  The  effects  of  reducing  clutter  will 
have  a  significant  effect  on  the  detection  of  targets  with 
small  radial  velocities. 

The  doppler  shift  in  a  bistatic  signal  return  is  not 
the  same  as  the  monostatic  return.  The  equations  are  more 
difficult  because  the  transmitter  and  receiver  are  not 
co-located.  The  bistatic  doppler  shift  is  a  result  of  a 
shift  due  to  the  radial  velocity  of  the  target  to  the  trans¬ 
mitter,  and  a  second  shift  due  to  the  radial  velocity  of  the 
target  with  respect  to  the  receiver. 

The  sol u  ions  to  the  bistatic  range  and  dc^pler  equa¬ 
tions  result  in  solution  families.  These  solution  families 
are  more  clearly  expressed  in  graplucai  form  using  Isorange 
and  Isodoppler  contours.  The  Isorange  contour  (Figure  2)  is 
the  family  of  possible  solutions  to  a  radar  range  return. 
Without  any  prior  knowledge  of  the  transmitter /receiver /tar¬ 
get  configuration,  the  isorange  contour  indicates  possible 
locations  for  the  target.  The  isodoppler  contours  (Figure  3) 
are  the  farily  of  constant  velocity  curves  (as  a  result  of 
doppler  shift) . 


Advantages 

Some  advantages  of  multistatic  systems  are  apparent 
when  applied  to  information  gathering  environments,  Advan- 
-• — ta9eS'4ncIu<i«t  -  -  -  - - - - - - — - - - -  . 
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1.  Undetectable  position  of  the  Receiver 

2.  Accuracy  in  Ranging 

3.  Operation  of  the  Receiver  in  Sanctuary 

4.  Cost  Savings 

An  airborne  bistatic  system  would  offer  some  surveil-' 
lance  advantages  in  that  the  radar  receiver  would  not  broad¬ 
cast  its  presence.  Transmitter  detection  of  the  receiver  is 
a  function  of  the  transmitter's  power,  coupled  with  the 
receiver's  distance  (two  way),  the  radar  cross  section  (RCS) 
of  the  receiver.  The  receiver  could  be  outside  tbe  trans¬ 
mitter's  detection  region  and  still  be  well  within  the 
exploitation  region. 

Farina  and  Hanle  have  demonstrated  that  a  multistatic 
ranging  system  is  more  accurate  provided  the  target  is  within 
the  triangulation  perimeter  (3*517).  The  airborne  bistatic 
radar  system  could  operate  in  the  confines  cf  a  netted  sys¬ 
tem  but  would  require  multiple  receivers  to  do  triangulation 
to  improve  accuracy.  In  general,  the  bistatic  radar  system 
will  not  be  as  accurate  as  a  monostatic  system. 

If  the  radar  could  operate  at  great  ranges  (100-200  KM), 
the  receiver  is  operating  in  a  relatively  safe  zone  or  sanc¬ 
tuary.  The  limitations  of  the  distance  between  the  trans¬ 
mitter  and  receiver  imposed  by  power,  coupled  with  the  lim¬ 
itations  imposed  by  waveform,  describe  the  safe  zone  oper¬ 
ating  limits. 

Whi  le  the  receiver  would  r^uir^  _ _ 

resources,  the  savings  in  design  costs,  by  eliminating  the 
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waveform  generator  portion  of  the  airborne  syst.em»  is  an 
attractive  trade-off  option.  Additionally*  the  savings  in 
receiver  platform  fuel  consumption  due  to  weight  elimination 
and  power  generation  could  be  attractive. 

Summary  of  Current  Research 

The  Bistatic  Radar  Problem  is  a  complex,  multifaceted 
arena  with  many  potential  areas  for  investigation.  Current 
issues  being  studied  by  the  technical  community  includei 
Bistatic  Cross  Sections,  Scintillation  Effects,  Polarization 
Effects  of  Returns,  Synchronizaticai  Techniques,  and  Signal 
Acquisition  Techniques. 

The  bistatic  radar  cross  section  is  not  the  same  as  the 
monostatic  radar  cross  section.  The  monostatic  signal 
return  is  a  function  of  the  target  distance  and  cross  sec¬ 
tion.  The  bistatic  cross  section  is  generally  larger  because 
of  the  larger  forward  energy  reflection  (15»3).  The  cross 
section  returns  are  dependent  upon  the  type  of  material  and 
the  type  of  background.  Many  technical  communities  are 
working  to  develop  models  to  analytically  predict  the  returns 
due  to  targets,  ground,  sea,  and  terrain.  Specific  studies 
representing  key  areas  of  interest  includei 

1.  RCS  Prediction  of  Low  Observable  Coated,  Objects  (5) 

2.  Ground  Clutter  (13) 

3.  Chaff  (8) 

4.  SeaAerrain  Scattering  Coefficients  (4) 
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A  Bistatic  Synthetic  Aperture  Radar  (SAR)  Self  Synchro¬ 
nization  Technique  has  been  proposed  by  Westinghouse  as  a 
result  of  their  studies  of  bistatic  radar  platforms  (1), 

The  report  outlines  a  synchronization  method  which  could  be 
used  between  a  receiver  and  a  cooperative  host.  The  concept 
of  a  Bistatic  SAR  was  proposed  and  demonstrated  in  the  T-Bird 
program  (Morenno) .  The  T-Bird  Program  demonstrated  the 
ability  to  gather  data  via  a  separate  airborne  tramsmitter 
receiver  pair. 


^ ^ ^ •  Pevelopment  of  the  Bistatic  Model 


Defining  the  Parameters 

The  bistatic  model  used  in  this  study  assumes  the  use 
of  a  pulse  doppler  radar  system.  The  model  assumes  the 
receiver  hardware  is  capable  of  detecting  and  analyzing  the 
radar  return.  The  parameters  used  in  the  model  are  the  com¬ 
mon  variables  affecting  any  pulse  doppler  signal  processing 
system.  The  model  consists  of  a  computer  program  which  uses 
the  parameters,  formulas,  and  concepts  explained  in  this 
chapter  (see  Appendix__A  for  a  complete  listing) .  The  com¬ 
puter  program  assumes  the  default  parameters  of  the  E3A 
AWACS  Waveform  but  can  be  modified  interactively.  The  sys¬ 
tem  parameters  include: 

Pulse  Repetition  Frequency  (PRF) 

Number  of  distinct  carrier  frequencies  used  (NPRF) 

Pulse  Duration  (TPULSE) 

Carrier  Frequency  (FO) 

Signal  to  Noise  Detection  Ratio  (SNR) 

Transmitter  Power  (P) 

Gain  of  the  Transmitter  Antenna  (GT) 

Gain  of  the  Receive  Antenna  (GR) 

System  losses  due  to  physical  construction  (LOSS) 
Receiver  Minimum  Signal  to  Noise  Detection  Ratio  (SNR) 
Average  Target  Radar  cross  section  (RCS) 

Receiver  Temperature  (T) 

Receiver  Bandwidth  (BW) 

Signal  to  Clutter  Power  Ratio  (SCR) 

Doppler  Frequency  Shift  of  the  Clutter  Cell  (FSC) 

Area  of  the  Clutter  Cell  (AC) 

Density  of  the  Clutter  Cell  (SO) 

Receiver  horizon  to  target  depression  angle  (RSEL) 
Transmitter  horizon  to  target  depression  angle  (XSEL) 
Receiver  baseline  to  target  azimuth  angle  (RSAZ) 
Transmitter  baseline  to  target  azimuth  angle  (XSAZ) 
Transmitter  Azimuth  Beam  width  (XSAZ) 

- Receiver  Azimuttr^eam-width -(RSAZ-) - -  - 

Transmitter  Elevation  Beam  Height  (XSEL) 


Receiver  Elevation  Beam  Height  (RSEL) 

I'lstatic  Angle:  Transmitter  to  Receiver  as  seen  from 
target  (A-B) 

Out  Of  Plane  Azimuthal  Angle  between  beams  (PHI) 

Receiver  Altitude  (RALT) 

Transmitter  Altitude  (XALT) 

The  bistatic  airborne  system  model  is  developed  in  the 
rectangular  coordinate  system  (X,  Y,  Z) .  The  Receiver/ 
Transmitter/Target  set  have  specific  physical  relationships 
(Figure  4).  The  receiver  antenna  is  the  base  for  all  meas¬ 
urements.  The  receiver  antenna  is  assigned  to  position 
X=0,  Y=0,  Z=  RALT  (RALT  default  is  9)  because  the 

receiver  is  assumed  to  be  airborne  at  a  constant  surveillance 
elevation.  The  target  can  be  anywhere  in  X,  Y,  Z  space. 

The  direct  path  between  the  transmitter  and  the  receiver  is 
known  as  the  baseline  in  the  bistatic  model.  Therefore,  the 
model  assumes  the  transmitter  is  located  at  a  point  along 
the  X  axis  (X  =  SEP)  and  at  an  altitude  Z  =  XALT.  This 
corresponds  to  position  SEP,  O,  XALT.  Note  that  all  dis¬ 
tances  used  in  this  model  are  in  kilometers. 

The  number  of  parameters  makes  the  assignment  of  mean¬ 
ingful  variables  difficult.  The  use  of  some  self  imposed 
rules  helps  to  aid  in  clarity.  The  transmitter  variables 
are  prefixed  with  X  and  the  target  variables  are  prefixed 
with  T.  For  example,  the  transmitter  azimuth  angle  is  coded 
XAZ  while  the  target  azimuth  angle  is  coded  TAZ. 

The  basis  for  the  mathematical  development  considers  a 
receiver  at  reference  position  0,  0,  RALT  (X  coordinate. 


Figure  4  Bistatic  Coordinate  System 

Y  coordinate,  Z  coordinate);  a  transmitter  at  position  SEP, 
0,  XALT;  and  a  target  anyvhere  in  XYZ  space  such  that; 

-50  <  X  -  250  kilometers 
-125  <  Y  <  125  kilometers 
0  <  Z  <  30  kilometers 

Radar  Range  Equation .  Classical  analysis  of  radar 
problems  begins  with  the  radar  range  equation  (12il9).  The 
equation  contains  the  parameters  of  interest  which  are 
affected  by  any  radar  system. 

^4 _ P  GT  GR  C^  RCS 

. . (WT^  T0^~BtrX“SNR~T -K - - 
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=  Power  of  the  Transmitter  (Watts) 

Gain  of  the  Transmitter  Antenna 

=  Gain  of  the  Receive  Antenna 

=  System  Losses  Due  to  Physical  Construction 

=  Receiver  Minimum  Signal  to  Noise  Ratio  Detection 

=  Speed  of  Light  (3  X  10®  Meters  per  Second) 

-23 

=  Boltzmann's  Constant  (1.38  X  10  Joules  per 
Degree  Kelvin) 

=  3.1416 

=  Average  Target  Radar  Gross  Section  (Square  Meters) 
=  Receiver  Temperature  (Degrees  Kelvin) 

=  Receiver  Bandwidth  (Hertz) 

=  Transmitted  Frequency  (Hertz) 

=  Distance  from  Transmitter/Receiver  to  Target 
(Meters) 


Note  that  the  radar  range  equation  can  be  rewritten: 


P  GT  C 


4  V  R1  FO 


4  Tr  R2 


4  TT  BW  L  T  SNR 


<  Transmitter  >  <  Target  >  <  Receiver  > 

Note  the  separation  of  the  parameters  into  transmitter,  tar¬ 
get,  and  receiver  dependent  areas.  These  are  the  parameters 
which  are  varied  to  examine  the  bistatic  case. 

Bistatic  Geometry.  The  geometry  of  the  bistatic  radar 
system  i-s  ellipsoid.  The  transmitter  and  the  receiver  are 
the  foci  of  the  ellipsoid.  The  surface  of  the  ellipsoid 
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forms  a  family  of  solutions  to  the  isorange  equation  prob¬ 
lem.  The  general  equation  for  an  ellipsoid  (13i51)  is*  , 

(  +  (  V|Vo  ^  ,  1  (3.2) 

In  this  case,  the  parameters  B  and  C  are  equal  because 
the  ellipsoid  has  no  eccentricity.  The  values  of  A  and  B 
are  developed  later  in  the  minimum  and  maximum  isorange  solu¬ 
tion  section. 

Recall  from  Figure  4  that  the  model  requires  the 
receiver  to  be  positioned  at  point  0,  0,  RALT  and  the 

transmitter  to  be  positioned  at  some  point  along  the  X  eixis 
specified  as  SEP,  0,  XALT.  If  the  transmitter  is  at  the 
same  altitude  as  the  receiver  (XALT  -  RALT),  equation  3.2 
becomes : 

(  ,%  (  v,^  )2  ,  J  (3.3, 

A  B  B 

The  detection  region  of  the  bistatic  system  is  limited 
by  the  power  of  the  transmitter.  The  transmitted  waveform 
must  have  enough  power  to  travel  (Figures  4  and  6)  a  dis¬ 
tance  from  the  transmitter  to  the  target  (Rl)  and  then  a 
distance  from  the  target  to  the  receiver  {R2)  and  still 
pro'/ide  a  large  enough  signal  to  noise  power  ratio  (SNR)  to 
insure  accurate  detection.  The  maximum  detection  range  is 
given  by  a  set  of  curves  referred  to  as  the  Ovals  of  Cassini 
as  shown  in  Figure  5.  Any  target  inside  the  oval  provides 
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Figure  5  Ovals  of  Cassini 


at  least  the  minimum  signal  to  noise  ratio  at'  the  receiver 
antenna.  Note  that  it  is  possible  to  have  two  completely 
disjoint  regions  of  detection.  The  transmitter  power  affects 
the  detection  region.  The  solution  set  to  the  ovals  of 
Cassini  will  yield  either  two  distinct  regions,  two  adjacent 
regions,  or  one  region. 

The  bistatic  radar  cross  section  of  a  target  varies 
from  the  monostatic  cross  section.  Investigation  into  the 
error  between  monostatic  and  bistatic  cross  section  has  shown 
that  the  Bistatic  Nonostatic  Equivalence  Theorem  is  valid 
only  for  small  bistatic  angles  (1S«3).  The  model  allows  for 
changing  the  RCS,  but  offers  ho  'monostatic  bistatic  compari- “ 


■T^he  bistatic  radar  range  equation  varies  from  the  mono¬ 
static  by  simple  substitution  of  Eq  (3.5)  into  Eq  (3.1). 

The  power  radiated  from  the  transmitter  spreads  in  a  sphere. 
The  gain  of  rhe  transmitter  antenna  accounts  for  the  direc¬ 
tionality  of  the  density.  The  radiated  power  has  a  specific 
density  which  is  a  function  of  the  square  of  the  distance 
from  the  target  to  the  transmitter.  At  that  point,  the  cross 
section  of  the  target  reflects  some  of  the  power.  The 
reflected  power  acts  as  if  the  target  were  an  isotropic 
radiator.  The  reflected  power  spreads  in  a  sphere.  The 
density  of  the  reflected  power  is  a  function  of  the  square 
of  the  distance  between  the  target  and  the  receiver. 


=  Rl^  R2^  (3.5) 

Oi serve  that  the  maximum  range  (R^)  is  a  constant.  The 

4 

Range  constant  (R  )  is  a  product  of  the  squares  of  the  trans¬ 
mitter  to  target  range  (Rl)  and  the  target  to  receiver  range 
(R2).  Therefore,  the  bistatic  radar  range  equation  becomesj 


7  2 

Rr  R2^  - 


P  GT  GR  C  RCS 
^4Tf)^  FO^  BW  L  SNR  T  K 


(3.6) 


The  transmitter  and  target  coordinates  and  tl«  Pythag¬ 
orean  Theorem  combine  to  show* 


Rl^  -  (SEP  -  X)^  +  +  (Z  -  XALT)^  (3.7) 
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Applying  the  same  to  the  target  receiver  pair* 


R2^  =  +  (z  -  raitV 


(3.8) 


Therefore* 

R'^  *=  [(SEP  -  X)^  +  y^  +  (Z  -  XALT)^] 


[x^  +  y^  +  (Z  -  RALT)^]  (3.9) 

Ability  of  System  to  Perf orm .  The  Bistatic  Radar  Iso¬ 
range  curves  (Figure  6)  are  calculated  by  extracting  informa¬ 
tion  from  the  radar  return  signal.  The  specific  measure¬ 
ments  are  the  Distance  (SEP)  between  the  transmitter  and  the 
receiver,  the  Angle  (ON)  between  the  transmitter  and  the 
target  as  measured  from  the  receiver,  and  the  time  difference 
(DT,  Figure  7)  between  receiving  a  direct  signal  (transmit¬ 
ter  to  receiver)  and  the  next  large  return  (transmitter  to 
target  to  receiver) . 

The  geometry  of  the  airborne  bistatic  radar  system 
yields  ellipsoidal  isorange  curves.  The  Hey  to  understand¬ 
ing  the  equations  is  to  have  a  consistent  frame  of  reference. 
Although  the  system  may  consist  of  three  airborne  entities 
(target,  receiver,  and  transmitter),  the  coordinate  system 
is  defined  by  making  the  receiver  the  reference  point.  This 
point  will  be  relative  in  XYZ  space  by  always  viewing  the 
receiver  as  being  at  the  origin  of  the  X  and  Y  axes  and  at 
_altitude.._Z  =_R^:T„KM._  {P,_Oj.RALT),.  All  distances  will  be 
measured  with  respect  to  the  receiver.  The  ’•BASELINE"  is 
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Fiqure  6  Elliose  of  Constant  Ranae 


Pulse  Train  at  Iranshitter  -^!  Uidth 


the  imaginary  line  between  the  receiver  and  the  transmitter. 
The  baseline  lies  along  the  positive  X  axis.  Therefore  the 
transmitter  is  always  located  at  some  separation  distance 
and  some  altitude  point  (SEP,0,XALT) . 

The  time  required  for  a  signal  to  travel  from  the  trans 
mitter  to  the  target  and  then  to  the  receiver  is  DT  seconds 
longer  than  the  time  required  to  travel  directly  from  the 
transmitter  to  the  receiver.  The  relationship  of  the  trian- 
gular  distances  isi 

R1  +  R2  =  C  *  DT  +  SEP  (3.10) 

This  result  is  a  family  of  solutions  (R1  +  R2)  which 
traces  out  an  ellipsoid.  The  point  on  the  surface  of  the 
ellipsoid  is  a  solution  of  interest.  This  family  of  solu¬ 
tions  is  called  an  isorange.  Obviously,  a  different  DT 
measurement  results  from  a  different  isorange  curve. 

Figure  7  points  out  several  areas  of  concern.  As  DT 
grows  small,  separdtion  of  the  pulses  (TV.IN)  in  the  pres¬ 
ence  of  noise  becomes  extremely  difficult.  Even  if  there 
were  an  absence  of  noise,  separation  of  the  pulses  is  dif¬ 
ficult  as  TWIN  approaches  zexo.  At  sotae  point,  the  two 
pulses  will  overlap.  The  receiver  is  unable  to  extract  any 
information  from  the  signal.  This  occurs  when  the  distance 
between  the  target  and  receiver  grows  small  (SEP  =  R2)  or 
when  the  bistatic  angle  (A-3)  approaches  160  degrees.  The 
receiver  must  be  able  to  detect  two  separate  signals.  This 


blind  area  is  a  disadvan^-.age  of  the  bistatic  model.  The 
blind  zone  is  treated  in  more  depth  in  the  ISOMIN  Develop¬ 
ment  Section. 

The  receiver  antenna  is  assirmed  to  be  a  directional 
scanning  antenna.  The  receiver  knovs  vhere  it  is  looking  at 
all  times.  Iherefcre,  when  a  signal  is  detected*  the 
bistatic  angle  (A-B)  is  known  within  the  limits  of  the  antenna 
beam  shape.  If  multiple  beams  are  used,  a  voliame  of  space 
can  be  scanned  with  greater  resolution.  Knowing  the  bistatic 
angle  {A-B)  gives  a  specific  solution  to  the  isorange  family. 
The  accuracy  of  solution  is  directly  dependent  upon  the  beam 
width  of  the  antenna  which  in  turn  limits  the  ability  to 
separate  the  target  and  transmitter  returns .  Several  alter¬ 
nate  solutions  have  been  offered  for  solving  this  problem 
(9:1) . 

One  possible  solu.ion  comes  from  basic  geometry.  This 
short  derivation  is  provided  to  show  one  possible  method 
for  receiver  calculation  to  target  location.  Knowing  (A-B) 
the  calculation  of  R1  and  R2  is  possible.  Therefore,  the 
target  range  and  direction  can  be  determined.  R2  is  derived 
from  the  Law  of  Cosines; 

R3  =  V  SEP^  +  (XALT  -  RALT)^ 

Rl“  =  R2^  +  R3^  -  2  R2  R3  CCS  (A-B 

and  because  the  bistatic  relationship  to  a  family  of  solu¬ 
tions  , 


R1  +  R2 

then, 

R1 

{K8  -  R2)^ 

K8^  -  2  K8  R2  +  R2^ 
Rearranging, 

R2  (2  R3  COS  (A-B)  -  2  K8) 
and  finally, 

R2 


=  Isorange  Constant  (K8) 

=  K8  -  R2 

=  R3^  >  R2^  -  2  R2  R3  COS  (A-B) 
=  R3^  +  R2^  -  2  R2  R3  COS  (A-B) 

2  2 
=  R3^  -  K8^ 

=  (R3^  -  K8^)  /  (2  R3  COS  (A-B) 
-2  K8) 


Limiting  Reqi ons 

liie  Rover  Limited  Region  ♦  Note  that  the  general  solu¬ 
tion  to  Eq  (3.9)  is  greatly  simplified  by  limiting  the  obser¬ 
vation  to  a  Sf^cific  set.  In  the  specific  case,  the  para- 
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meters  SEP,  Z,  XALT,  and  R  are  constants  which  are  depend¬ 
ent  on  the  transmitter  receiver  target  geometry  and  hardware. 
Equation  (3.9)  is  solved  for  a  specific  Y  in  terras  of  X  by 
use  of  the  quadratic  formula  (13:32).  The  solution  may 
yield  one  of  the  three  region  shapes  shown  in  Figure  5 
depending  upon  the  values  assigned  to  the  constants . 

Let  Cl  =  (Z  -  RALT)^  (3.11) 


C2  =  (Z  -  XALT)^ 


(3.12) 


Then,  from  (3.9): 


r'^  =  [(SEP  -  X)^  +  +  C2]  [X^  +  +  Cl] 


(3.13) 


Expanding  the  terms  and  rearranging  yields: 

0  =  +  (X^  +  (SEP-X)^  +  C2  +  Cl)  -P"  +  X^ 


((SEP-X)^  +  Cl))  +  Cl  [(SEP-X)^  +  C2]  -  r"^  (3.14) 
Note  that  this  is  in  quadratic  form: 


AY"^  +  BY^  +  C  =  0 


r2  _  -  B±  »Vb"^-4AC) 


And  therefore: 


Y  =  ±  •  V  Y^ 


Where: 


A  =  1 


B  =  +  (SEP-X)^  +  C2  +  Cl 


C  =  X^  [(SEP-X)^  +  C2]  +  Cl  C(SEP-X)^  +  C2]  - 


C3  =  C2  +  Cl  +  SEP 


(3.15) 


CT  ="  -""2  SEP 


.  (3:16) 


.^v 


G5  >  Cl.^C4 


^^3vl7) 


^ci  -  s#^_  +  ;^ci;?c2^ 


^(3.i8) 


B4  =  -  X  SEP  +03/2 


:<3,19) 


CC  =  X^  +  C4  X^  +  C3  X^  +  C5  X  +  G6  ;  ^(3.20) 


Then  solving  for  solutions  of  Yj 


Y  =  ±  -  V  -  B4  +  V  {-B4^  -  C3) 


(3.21) 


Isorange  Formula  Development .  It  was  shown  earlier 
that  the  solution  set  to  the  bistatic  range  equation  resulted 
in  a  family  of  solutions.  The  family  of  solutions  occupy  an 
ellipsoidal  space.  That  is  to  say,  the  solution  set  is  the 
surface  of  an  ellipsoid.  Recall  Eq  (3.2).  This  equation 
can  be  solved  for  the  roots  of  one  variable  (X,  Y  or  Z). 

For  this  model,  a  plane  of  interest  is  specified  (Z  =  Z*). 
The  general  solution  set  for  Y  in  terms  of  X  isi 


,  I  .  1  .  (  5.^1^  .  (  X..-  SJP/2  ,2  ,3  ^2, 

B  BA 


Then, 


-  (2  -  RALT)^  -  (  °  (3.23) 

A 


or,  considering  both  roots  as  solutions i 


B  (X  -  SfiP/2) 


This  general  equation  holds  only  for  an  ;airrporhe  to^ 
airborne  transmitter  receiver  pair  at  the  sarae^  a.ititude^\^^^ 
is  more  realistic  to  assume  that  the  non*“cooj^rttiv6  t 
mitter  could  be  at  an  altitude  of  up  to-  30  kilometers,  .  ; 

Allowing  for  this  possibility  adds  complexity  to  the  eg^- 
tion  but  only  through  a  rotation  and  translation -of  the.k  ‘ 
and  Z  axes. 
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Figure  8  Rotation  of  the  Bistatic  Model 

The  rotation  angle  is  THETA  =  ARCSIN  (XALT  -  HALT)  / 
SEP.  This  makes  the  coordinates  of  the  new  axis  X  =  RALT 
SIN  THETA,  Y  =  0,  Z  =  RALT  (1  -  COS  THETA).  The  trans¬ 
forming  equations  (13 >36)  are  then: 


X  =  X'  cos  THETA  -  Z 'SIN  THETA 


+  RALT  SIN  THETA 


Z  =  X’SIN  THETA  +2 'COS  THETA 


0.25)  r 


m 


+  RALT  (1  -  COS  THETA) 


(3,26) 


Substitution  of  these  transforms  into  the  general  equation 
yields: 


-  (X‘  SIN  THETA  +  Z'COS  THETA  +  RALT  COS  THETA)^ 

-  (B/A)^  (X'COS  THETA  -  SIN  THETA  (Z  -  RALT) 

-  SEP/2)^  (3.27) 


Assigning  the  following  simplifying  variables  to  the  con¬ 
stants  and  dropping  the  prime  notation  gives: 


(3.28) 


SIN  THETA  =  STHETA  =  (XALT  -  RALT)  /  SEP  (3.29) 

COS  THETA  =  CTHETA  =  COS  (ARCSIN  THETA)  (3.30) 

CIO  e  Z  COS  THETA  -  RALT  COS  THETA  (3.31) 

Cll  =  Z  SIN  THETA  -  RALT  SIN  THETA  +  SEP/2  (3,32) 

The  equation  in  terms  of  X  and  Y  yields, 

-Y  ±  ,  --;/-b2  (x-stHETA~+“C10)^  '-  ”C12  ’TX-TTHETA"  -  ‘clll^^ - 


i 

; 


ft 


I 


'i 


i 


I 


;  Ispr^q e  B1  Ind  Zone,  ( ISQMIN )  »■  The .  minimxim  ' 

the- s  Glut  -trsuid-^ 

along  this  patl}  '(4j  5»  15)  .  Therefore,  a  starting  point 
must  be  found  at  .  some  practical  rartge .  logical  starting 

point  is  the  edge  of  the  “Isorange  Blind  Zonii*'. 

The  blind  zone  (Figure  9)  is  the  area  which  results 
from  the  receiver  being  unable  to  determine  where  the  pulse 
traveling  along  the  direct  path  between  the  transmitter  and 
the  receiver  stops  and  the  pulse  traveling  along  a  reflected 
path  starts.  This  blind  zone  will  be  a  function  of  the  pulse 
duration  (TPULSE).  The  minimum  blind  zone  isorange  is  twice 
the  pulse  duration  times  the  speed  of  light,  plus  the  separa¬ 
tion  distance  between  the  transmitter  and  receiver. 


R3  =  V  SEP  +  (XALT  -  RALT) 


(3.34) 


ISOMIN  =  (2  *  TPULSE  *  C)  +  R3 


(3.35) 


Beyond  the  bl ' nd  zone  the  isoranges  need  to  be  incre¬ 
mented  at  some  regular  interval  up  to  some  limiting  factor. 
The  model  limit  is  300  Kilometers.  Often  the  Pulse  Repeti¬ 
tion  Frequency  (PRF)  will  dictate  a  maximum  isorange. 

Maximum  Isoranqe  ( ISOMAX ) .  The  PRF  limitation  occurs 
when  the  transmitted  PRF  becomes  so  high  that  the  bistatic- 
J>'ilse_is  not  received  until,  during  or  after 
a  second  pulse  is  received  on  the  direct  path.  At  a  low  PRF, 
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Figure  9  Isorarxge  Detection  Region 

the  direct  pulse  is  received  folloved  some  time  later  by  the 
target  reflected  pulse.  However,  when  the  PRF  becomes  suf¬ 
ficiently  high,  the  time  of  reception  of  the  bistatic  pulse 
can  be  delayed  long  enough  so  that  it  overlaps  the  next 
direct  pulse.  The  overlap  begins  to  occur  when  the  leading 
edge  of  the  second  transmitted  pulse  arrives  before  the 
trailing  edge  of  the  first  bistatically  reflected  pulse  (4i9)i 

R1  +  R2  +  TPULSE  C  >  (C  /  PRF)  +  R3  (3.36) 

This  leads  to  the  exclusion  region  (Figure  9)  of  the  ellipse 
defined  byt 


R1  +  R2  =  (C  /  PRF)  +  R3  -  TPULSE  C  (3,^37) 

But,  if  multiple  PRFs  (NPRF)  are  used  by  the  transmitter, 
the  exclusion  region  can  extend  beyond  these  bounds  (I2ill4)* 
The  maximum  isorange  (ISOMAX  =  R1  +  R2)  is  therefore 
extended  toj 

ISOMAX  =  (NPRF  *  C  /  PRF)  +  R3  -  TPULSE  C  (3.38) 

Incrementing  the  Isoranges.  It  is  often  necessary  to 
draw  several  isorange  curves  between  the  minimum  (ISOMIN) 
and  the  maximum  (ISOMAX)  isorange  curves.  Choosing  10 
Kilometers  as  an  interval,  the  first  isorange  is  needed  such 
that  ISOMIN  <  ISO(I)  <  ISOMAX. 

Letting, 

Q  =  INT  [  ISOMIN  /  10  ]|  Greatest  Integer  (3.39) 

ISO(I)  =  I  *  Q  *  10  (3.40 

I  is  the  incremented  from  1  to  some  number  (K)  such  that 
ISO(I)  does  not  exceed  ISOMAX  or  300  Kilometers  (limit  of 
the  model).  For  each  ISO(I),  a  unique  solution  set  for  the 
equation  is  determined.  Frcwi  the  properties  of  the  ellipsoid 
(I3j38)j 


A  =  IS0(I)/2 


(3.41) 


I 
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and  because 


'■  .'  tv' 


;■■;■  ,.■■  ..  .,  ’•r  '•■;:■;* ■.■:'•\-..•■■^:.:■VV,■>T^^ 


B  «  V  -  ( SEP/2 


(3.43) 


term  A  simplifies  to: 


A  =  (NPRF  *  C  /  PRF  +  R3  -  TPULSE  C)  /  2  (3.44) 


Likewise,  B  reduces  to: 


B  =  V  [[(NPRF  *  C  /  PRF)  -  TPULSE  C]  *  [(NPRF  *  C  /  PRF) 


+  (2  *  R3)  -  TPULSE  c]]  /  2 


(3.45 


Note  that  in  most  applications,  TPULSE  «  1/PRF.  This 


simplification  holds  for  most  systems  as  the  prf  is  reduced 


to  extend  the  system  visibility.  The  zones  are  then  sepa~ 
rated  by  ellipses  whose  characteristics  A,B  reduce  to: 


A  =  ((NPRF  *  C  /  PRF)  +  R3)  /  2 


(3.46 


B  *=  V  CCnPRF  *  C  /  prf]  *  [{NPRF  *  C  /  PRF)  f  (2  *  SEP)]]  /  2 


For  example,  to  determine  A  and  B  for  ISO(I)  =  ISOMAX,  the  ■ 
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The  doppler  shift  of  a  signal  for  the  mpn^tatic  '4  -:  ^ 

tern  is  a  function-  of-''^he'ifriMVJBncy'-of  ' trensnu.sipipn,-’aha^;?J^p'ffE^%j«'%^%ea5 

■  '  "  ;  r 

velocities  of  the  transmitter,  target,  and  receiver.  The 


velocity  of  the  target  is  alvays  meastired  radially  fr(Mn  the 
transmitter  to  the  target.  The  radial  velocity  causes  a 
positive  (target  and  receiver  are  closing),  negative  (target 
and  receiver  are  separating)  or  zero  (target  and  receiver 
remain  at  a  constant  distance)  shift  in  frequency. 

The  classical  doppler  model  is  based  on  a  single  frame 
of  reference  (transmitter  and  receiver  are  co-located) .  The 
doppler  shift  of  a  transmitted  vave  propagating  through  space 
from  the  transmitter,  striking  an  object  (the  target),  and 
being  returned  to  the  receiver  is  twice  the  radial  velocity 
of  the  target  divided  by  the  transmitted  vave  length.  The 
shift  in  the  frequency  between  the  transmitted  signal  and 
the  returned  signal  is  the  doppler  shift  (12i69).  Analytic- 
allyt 

FD  =  2  V  /  L  (3.48) 


Where t 

FD  =  Doppler  Shift  (Hertz) 

V  =  Radial  Velocity  Between  Tvo  Objects  (Meters  per 
Second) 


L  =  Wave  Length  of  the  Transmitted  Frequency  (Meters) 
The  frequency,  vave  length,  and  speed  of  light  relationship 


i 
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C  =  L  FO 


(3.49) 


Therefore,  . 


FD  -  2  V  FO  /  C 


(3.50) 


The  bistatic  doppler  model  is  not  as  easily  derived* 
This  results  because  the  transmitter  and  the  receiver  are 
not  located  at  the  same  point  in  space.  In  fact,  it  is  con¬ 
ceivable  to  imagine  the  transmitter  receiver  pair  separated 
by  large  distances  (several  hundred  kilometers).  The  net 
doppler  shift  of  the  transmitted  frequency  (FF)  is  a  func¬ 
tion  of  the  radial  velocity  of  the  transmitter  with  respect 
to  the  target  and  a  second  doppler  shift  due  to  the  radial 
velocity  of  the  target  with  respect  to  the  receiver. 

The  transmitter  will  broadcast  a  signal  at  a  given 
Frequency  (FO).  The  target  will  reflect  the  signal  at  same 
doppler  shifted  frequency  due  to  the  radial  velocity  of  the 
transmitter  target  pair  (FDl).  The  reflected  signal  can  be 
modeled  as  if  it  were  an  isotropically  radiated  signal  at 
scane  new  frequency  (FN  =  FO  FDl).  The  signal  will  be 
acquired  at  the  receiver  at  some  second  doppler  shifted  fre¬ 
quency  (FF)  which  is  a  function  of  the  radial  velocity  of 
the  target/receiver  pair  (VI). 


FD  =  VI  /  L 


(3.51) 


=  VI  FO  /  C 


(3.52) 
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In  the  model,  the  transmitter  and  target  velocities  are 
specified  as  a  vector  with  magnitude  and  direction.  The 
velocities  are  user  specified  as  a  magnitude  (VT,VX), 
some  angle  of  azimuth  heading  (TAZ,XAZ)  and  some  angle 
depicting  rate  of  climb  (TEL,XEL).  The  azimuth  direction 
is  always  referenced  to  the  baseline  which  is  the  direct 
path  between  the  transmitter  and  receiver.  The  receiver  is 
assumed  to  be  flying  a  constant  elevation  surveillance  orbit. 
The  transmitter  will  also  have  a  zero  angle  for  TEL  (but  is 
not  model  limited) .  The  transmitter,  whether  ground  based 
or  airborne,  is  assumed  to  be  very  limited  in  maneuverability 
due  to  the  physical  characteristics  (construction)  of  the 
antenna.  An  obvious  exception  would  be  a  Phased  Array  Radar 
which  will  not  be  considered  he.e. 

The  radial  velocities  of  the  transmitter/target  (target/ 
receiver)  pair  are  computed  for  any  X,Y,Z  point  using  vec¬ 
tors.  The  radial  velocity  is  the  vector  dot  product.  The 
velocity  vector  is  dotted  into  the  unit  vector  separating 
the  two  objects. 

V  =  A  •  B 

The  velocity  vectors  VI,  V2  are  decomposed  into  X,Y,2 
components  for  ease  of  operation. 

\OCZ  =  VX  Sin  XEL  ;  Z  component  of  Transmitter  (3.53) 


\ork'  =  VX  Cos  XEL  S in  XAZ  j  V  component  of  'Transmit ler 


VXX  =  VX  Cos  XEL  Cos  XA2  ;  X  component  of  Transmitter 

(3.55) 

VTZ  =  VT  Sin  TEL  j  Z  component  of  Target  (3,56) 
VTY  =  VT  Cos  TEL  Sin  TAZ  ;  Y  component  of  Target 

(3.57) 

VTX  =  VT  Cos  TEL  Cos  TAZ  ;  X  component  of  Target 

(3.58) 

The  general  form  of  the  dot  product  (16j635)  is: 

A  ‘  B  =  ai  bl  +  a2  b2  +  a3  b3  (3.59) 

The  vector  from  the  receiver  to  the  general  point  XYZ  is: 

A  =  X  i  Y  j  +  (Z-RALT)  K  (3.60) 

The  vector  from  general  point  Xi'Z  to  the  transmitter  is: 

B  =  (SEP-X)  i  -  Y  j  +  (XALT-2)  k  (3.61) 

Both  v'^ectors  A  and  B  must  be  normalized.  The  normaliza¬ 
tion  is  reguired  as  only  the  unit  position  vector  is  required 
Therefore: 

Let  A1  ^  ♦  (Z-RALT)^  (3.62) 

and 

Let  B1  =  y!  (SEP-X)“  +  Y^  (XALT-Z)'^  (3.63) 


Recall  that  the  vector  between  two  points  has  direction 
as  veil  as  magnitude.  Because  the  vector  from  X,  Y,  Z  to 
SEP,  0,  XALT  is  the  negative  of  the  vector  from  SEP,  0,  XALT 
to  X,  Y,  Z,  a  factor  of  -1  must  be  included  in  the  trans¬ 
mitter  component  to  insure  the  doppler  shift  is  negative 
vhen  separating  and  positive  when  closing.  The  radial  veloc¬ 
ity  of  transmitter  with  respect  to  the  target  due  to  the 
velocity  of  the  transmitter  alone  is  thens 

VIR  =  [(X  -  SEP)  VXX  +  Y  VXY  +  (Z  -  XALT)  VXZ]  /  B1 

(3.64) 

The  radial  velocity  of  the  transmitter  with  respect  to 
the  target  due  to  the  velocity  of  the  target  alone  is  then* 

V2R  =  [(SEP  -  X)  YTX  -  Y  VTY  +  (XALT-Z)  VTZ'J  /  B1 

(3.65) 

The  total  radial  velocity  due  to  the  transmitter  target 
pair  is*. 

VTQT  =  VIR  V2R  (3.66) 

The  doppler  shift  is  a  function  of  the  frequency  and  wave 
leng’cht 

V  =  L  FD  (3.67) 

FDl  =  VTOT  /  L 

or,  by  combining  and  simplification: 


Lai; 
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FDl  =  L(SEP-X)  (VTX  -  VXX)  +  Y  (VXY  -  VTY) 


+  (XALT-Z)  (VTZ  -  VXZ)  ]  /  (B1  L)  (3.68) 


At  this  point  the  new  radiated  frequency  is  the  transmitted 
frequency  plus  the  doppler  shift.  The  new  frequency  ist 


FN  =  FO  +  FDl 


(3.69) 


And  of  course,  this  new  frequency  has  a  new  wave  lengtht 


LN  =  C  /  FM 


The  radial  velocity  of  the  target  v’ith  respect  to  the 


(3.70) 


receiver  due  to  the  velocity  of  the  target  (model  assumes 
receiver  is  zero  reference)  is  then} 


V2R  =  [VTX  X  +  VTY  Y  +  VTZ  (Z-RALT)]  /  A1  (3.71) 


Using  the  same  development,  the  second  doppler  shift  is; 


FD2  =  V2R  /  LN 


=  (-1)  [VTX  X  +  VTY  Y  +  VTZ  (Z-RALT)]  /  (A1  *  LN) 


(3.72) 


The  net  doppler  shift  is; 


FF  =  [(SEP-X)  (v'TX  -  \OOC)  +  Y  (\DCY  -  VTY) 


(XALT-Z)  (VTZ  -  VXZ)'i  /  (B1  L) 


-  [VTX  X  +  VTY  Y  +  VTZ  (Z-RALT)]  /  (A1  LN )  (3.73) 


This  is  the  general  form  of  the  solution  space.  A  map 
of  the  isodoppler  contours  is  desired.  An  array  (TABLE)  is 
set  up  to  hold  the  general  solution  set  for  the  region  of 
interest.  A  plane  of  interest  is  examined  by  specifying 
the  Z  plane  (altitude)  desired.  For  example,  if  the  plane 
of  interest  is  the  doppler  contours  of  the  ground  plane 
(Z  =  0),  the  array  is  scanned  for  contour  information  and 
a  table  of  values  is  generated.  A  contour  map  (isodoppler 
lines)  is  constructed  by  a  hand  drawing  the  contour  lines 
on  the  map.  The  doppler  value  equations  are  nonlinear  and 
do  not  yield  to  a  simple  plotting  routine. 

Clutter 

Clutter  Parameters .  Clutter  is  the  unwanted  radar 
echo  from  an  object  or  background  (12:470).  Clutter  is 
usually  the  radar  return  from  a  ground  cell,  insects,  birds, 
or  possibly  a  storm  cell.  The  clutter  power  is  often  larger 
than  any  possible  target  power  signal .  This  occurs  because 
the  radar  cross  section  of  a  clutter  cell  is  much  larger 
than  that  of  a  probable  target.  The  radar  designer  is  con~ 
cerned  with  the  characteristics  of  the  unwanted  signal 
because  techniques  are  available  to  enhance  the  detection 
of  targets  if  the  clutter  can  be  removed.  Clutter  signal 
properties  can  be  predicted.  Therefore,  the  designer  builds 
"filters'*  to  modify  the  signal  returns  in  order  to  find  prob 
able  targets.  The  largest  clutter  return  comes  from  ground 


The  target  will  reflect  an  amount  of  power  proportional  to 
its  radar  cross  section.  The  reflected  power  has  an  intent 
sity  and  doppler  shift.  The  clutter  cell  "behind  the  target 
also  reflects  energy.  The  signal  power  (in  terms  already 
defined)  isj 


P  GR  GT  RCS 
(4it)^  R2^  FO^ 


(3.74) 


The  clutter  power  has  the  same  equation  with  slight  modifies 
tion.  The  surface  clutter  cross  section  is  defined  ast 

PCS  =  SO  AC  (3.75) 


Where: 

SO  =  Surface  Cross  Section  per  Unit  Area 

AC  =  Radar  Cross  Sectional  Area  of  Illumination 

The  bistatic  case  modifies  AC  to  mean  the  common  area 
between  the  transmitter  beam  and  the  receiver  beam.  Figure 
10  illustrates  the  common  area  for  several  possible  cells. 
The  Surface  Cross  Section  is  determined  from  a  model  which 
results  from  experimental  data  (15i20-23).  SO  varies  depend 
ing  upon  the  composition  of  the  surface.  For  exaitple,  sea 
and  forest  clutter  have  different  values.  Several  bistatic 
models  have  been  developed  to  explore  the  SO  phenomenon. 
References  (5)  and  (I5)  catalogue  many  models  and  much  of 
the  research  in  this  area . 


The  clutter  power  ist 


P  GR  GT  AC  SO  C?  ,(3,76); 

{4Tr)^  Rl^  R2^  FO^ 

Dividing  Eg  (3.78)  by  Eg  (3.80),  the  Signal  to  Clutter Jtatip 
is: 

SCR  =  - —  (3.77) 

SO  AC 

Clutter  Cell  Area.  The  clutter  area  is  dependent 
upon  the  beam  shapes  and  the  processing  eguipment.  The 
clutter  cell  is  either  beam  width  limited,  resolution  time 
limited,  or  doppler  freguency  limited  (15:12).  This  model 
uses  the  Resolution  time  limited  cell  area  developed  by 
Weiner  and  Kaplan  (15} 16).  The  bistatic  models  for  clutter 
assume  the  flat  earth  model.  The  clutter  cell  area  for  low 
grazing  angle  conditions  (RSEL  <  30  degrees)  is, defined  as: 


AC  = 


R1  XSAZ  SEC  (  A  -  B  ) 
2  BW  2 


;  XSAZ  <  RSAZ  R2/R1  «  2  RAD 


or; 


(3.78) 


AC 


^  R2_RS^SEC  (  A_-_B  )  ^  ^  Rl/^2  «  2  RAD 


2  BW 


Where: 

RSEL  =  Receiver  Horizon  to  Target  Depression  Angle 
(Degrees) 

XSAZ  =  Transmitter  Azimuth  Beam  Width  (Radians) 


RSAZ  =  Receiver  i^imuth  Beam  Width  (Radians) 


A-B  =  Angle  [Transmitter  to '  Recover  from 

Target]  (Radians) 

Using  the  Law  of  Cosines  and  trigonometric  identities,  it 
can  be  shown  that  the  termi 


SEC^  - 2 

^  2  R1  R2  +  Rl^  +  R2^  -  SEP  -  (RALT  -  XALT) 

(3.79) 

Therefore,  AC  can  be  estimated  in  terms  of  system  position. 

It  is  assumed  that  the  transmitter  beam  width  (XSAZ)  is 
known  or  at  least  a  good  estimate  can  be  made. 

Surface  Cross  Section  Per  Unit  Area,  Figure  11 
shows  an  estimate  of  the  Surface  Cross  Section  Density  Model 
(13:21),  The  graph  values  (for  3g  Hz  only)  have  been  modeled 
by  the  second  order  equation: 


SO  =  V  (X-90)  3.21111  -32  j  90  <  X  <  270 

(3.80) 

or: 


SO  =  V  (90-X)  10.000  -32  ;  0  <  X  <  90 

i  270  <  X  <  360 
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Weiner  has  shown  that  this  general  model  can  be  extrapolated 
to  other  grazing  angles  by  the  addition  of  a  correction  term 
(13j20)  to  the  value  obtained  from  the  graph.  The  correc¬ 
tion  term  ist 


COR  = 


PHI  ADJ 
180 


(3.81) 


Where* 


ADJ  =  -8  dB  I  (XSEL  +  RSEL)/2  <1.5  (3.82) 


or* 


ADJ  = 


9.5  Log,„  ((XSEL  +  RSEL)/2) 


(XSEL  +  RSEL)/2  >1.5 


Where j 


PHI  =  Out  of  Plane  Azimuthal  Angle  Between  Beams 
XSEL  =  Transmitter  Depression  Angle 
RSEL  =  Receiver  Depression  Angle 

Doppler  Spread  of  Clutter  Cell*  The  clutter  cell  power 
is  concentrated  in  the  frequency  domain  about  the  cell  dop- 
pler  shift  (Figure  12).  By  estimating  the  width  of  the 
clutter  cell  frequency  shift,  the  receiver  can  dynamically 
’’filter”  the  returned  signal.  The  ground  clutter  cell  dop- 
pler  frequency  shift  varies  as  the  geometry  of  the  common 
beam  area  changes.  The  equations  in  the  Clutter  section  are 
used  to  calculate  the  doppler  shift  of  any  specific  point. 

A  method  for  estimating  the  difference  in  the  maximum  and 
minimum  doppler  shifts  of  any  clutter  cell  is  required. 


FQ 


IN  BEAN  ClilirER  BOPFLER  SPREAD 


FREQUEHCy  DOMAIN 


Flgure~l‘2  cltrtter~Pover'  in  -the-Freguency-Domain 


The  model  assumes  the  use  of  scanning  and  illumination 
beams  vith  width  (degrees  in  azimuth)  and  height  (degrees  in 
elevation).  The  parameters  XWAZ,  XWEL,  RWAZ,  and  RWEL  , 
define  the  transmitter  and  receiver  beam  widths  and  heights 
respectively.  A  simple  method  for  determining  the  doppler 
spread  is  to  estimate  the  spread  based  on  the  beam  widths. 
This  is  done  by  calculating  four  points  for  each  beam.  These 
points  are  approximately  on  the  3dB  contour  of  the  beam 
(Figure  13).  A  distance  measure  is  made  to  each  of  the 
eight  points.  The  four  smallest  distance  measures  dictate 
the  four  points  used  to  estimate  the  doppler  frequency 
spread.  This  is  done  by  calculating  the  doppler  shift  to 
the  four  points  using  Eq  (3.77).  The  minimum  value  is  sub¬ 
tracted  from  the  maximum  value  to  obtain  the  estimate  of 
the  frequency  shift.  This  method  is  used  in  Subroutine 
DOPSD,  Appendix  A. 


Figure  13  Geometry  of  Common  Beam  Area 


Sununary 


The  equations  developed  in  this  chapter  are  used  in  the  ; 
Bistatic  Radar  Computer  Model.  In  this  chaptfEu:  the  ;iimits  ; 
of  the  model  were  defined.  These  limits  included  the  X,  Y,  Z 
region  of  interest  as  well  as  the  transmitter,  receiver,  and  . 
target  parameters.  These  equations  are  the  foundations  upon  , 
which  the  computer  model  was  built.  The  computer  model 
allows  the  operator  to  change  the  Vcdues  of  the  system  frcmi 
the  default  values.  Figure  14  shows  a  sample  menu.  Based 
on  the  values  entered  by  the  operator,  a  plot  is  generated 
(Figure  15).  This  plot  shows  the  isorange  and  power  limita- 
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Figure  14  Computer  Model  Menu 
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IV.  Exercising  the  Model 


Introduction 

This  chapter  outlines  a  generic  form  of  the  bistatic 
radar  model  and  demonstrates  the  effects  of  changing  the 
parameters  of  interest.  The  interactive  software  provides 
many  of  the  parameters  which  may  be  changed  by  the  user  when 
exercising  the  model.  This  chapter  demonstrates  the  effects 
of  transmitter  and  receiver  hardware  (power,  waveform,  etc.), 
platform  motion,  and  platform  orientation  on  the  bistatic 
detection  regions,  the  doppler  effects,  the  signal  to  clutter 
ratio  and  the  doppler  spread  of  the  received  signal.  Finally 
a  typical  scenario  is  illustrated.  For  all  cases  used  in 
this  section  a  generic  bistatic  set  is  assumed  and  changes 
are  shovn  relative  to  the  generic  set.  The  generic  model 
has  parameters  as  shown  in  the  sample  computer  menu  (Figure 
16).  The  parameter  units  and  range  of  acceptable  values  are 
displayed  when  changing  any  of  the  26  parameters. 

Limiting  Regions 

The  Power  Limited  Region.  The  detection  region  whicdi 
is  limited  by  the  transmitter  and  receiver  hardware  is 
described  by  the  Ovals  of  Cassini.  Figure  17  show^  the 
bistatic  detection  region  for  the  ground  plane  (2  **  0) . 
Decreasing  transmitter  power  causes  a  collapse  of  the  ovals 
as  demonstrated  in  Figure  18.  Finally,  Figure  19  shows  the 

^  . 
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igure  17  Gi^eric  Detection  Region  {Power  =  0*01  Megawatts) 


Figure  13  Detection  Region  (Power  =  0.0095  Megawatts) 


pout:r  ilimited  pegio 


Figure  19  Detection  Region  (Power  =  0.0090  Megawatts) 


detection  region  when  an  even  smaller  amount  of  power  is 
used  and  the  ov^is  divide  into  two  specific  regions. 


w 


The  shape  of  the  detection  region  changes  as  a  function 
of  the  platform  separation  distance  (SEP).  Figures  20  and 
21  display  the  change  in  the  detection  region  as  the  separa¬ 
tion  region  between  the  tran'^mi  tt-pr  and  ronc-^frar-  . 

150 -T— frlOT^TIiC  bCTCOTlOH  HLQlOttO - 1 - ! - 1 


POVJER  LIMITEDI  REGIO 


/  M 


rj.yure  zu  uerecrion  Kegion  i^eparation  =  190  KM) 


POKTR  LI.MSTTED  REGION 


Figure  2]  Detection  Region  (Separation  =  160  KM) 


Changing  any  of  the  parameters  which  effect  the  range 
constant  (r"^)  yields  similar  results.  Recall  from  Eq  (3.5) 
(3,6)  and  (3.9)  that  the  range  constant  is  a  function  of  the 
transmitter  and  receiv'er  hardware  as  well  as  the  target 
location  and  cross  section.  The  model  takes  into  considera¬ 
tion  varying  receiver  and  transmitter  altitudes.  Altitude 
changes  affect  the  geometry  of  the  model  by  tilting  the 
ovals  and  ellipsoids  producing  an  asymmetric  figure.  An 
example  of  an  airborne  to  ground  pair  is  shown  in  Figure  22, 
The  detection  region  varys  but  is  not  quite  as  dramatic  as 
the  change  in  the  waveform  limited  characteristic  regions 
defined  by  the  isorange  curves. 


POV.’ER  LIMITED  REGIO! 
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Figure  22  Airborne  Receiver  to  Ground  Transmitter 


Wavef orm  Limit ed  Regions . 

Eva.l via t ion  of  the  Blind  Zone  ( ISOMIK )  .  The  blind 
zone  is  the  region  between -the  transmitter  and  receiver* 


defined  by  an  ellipsoid,  where  tlie  receiver  cannot  detect 
any  target.  Heoall  from  the  ISOKIN  Development  Section 
Eq  (3.38)  that,  the  blind  zone  is  defined  by  the  isorange 
curve  (ISOMIN): 

ISOMIN  -  2*TPULSE*C+SEP 


The  ISOMIN  curves  shown  in  Figures  16  thru  22  are  for 
pulse  widths  of  5  microseconds.  It  is  easily  demonstrated 
that  a  longer  pulse  significantly  increases  the  width  of 
the  blind  zone.  Figures  17  and  23  demonstrate  the  effect  of 
increasing  the  pulse  width  from  5us  to  30us. 


Figure  23  Detection  Region  for  Pulse  Width  30  Microseconds 


The  width  of  the  blind  zone  ellipsoid  can  be  shown  to  be: 


WIDTH  =  SORT  (ISOMIN^  -  R3^) 


(4.1 


Table  I  lists  the  ISOMINs  for  a 

constant 

separation  . 

of  200  KM. 

TABLE  I 

Comparison  of  Pulse 

Width, 

Pulse  Width 

Isomin  and  Ellipsoid 

ISOMIN 

Width 

Ellipsoid  Widt 

(microseconds) 

(kilometers) 

(kilometers) 

.1 

200.06 

4.6 

1.0 

200.60 

15.5 

5.0 

203.00 

34.8 

10.0 

206.00 

49.4  .. 

15.0 

209.00 

60.6 

20.0 

212.00 

70.3 

25.0 

215.00 

79.9 

30.0 

218.00 

86.8 

When  the  data  is  put  into  graphical  form,  it  is  easy  to 
recognize  that  a  small  pulse  width  is  desirable  to  minimize 
the  blind  zone. 

Evaluation  of  the  Maximum  Range  (ISOMAX) .  The 
maximum  visibility  of  the  bistatic  system  was  shown  to  be  a 
function  of  the  pulse  repetition  frequency  (PRF),  the  nxamber 
of  distinct  pulse  repet ion  frequencies  (NPRF),  and  the  pulse 
width  (TPULSE),  Examination  of  Eq  (3,44)  yields  the  follow¬ 
ing  observations:  increasing  the  pulse  width  will  lower  the 
maximum  range  (compare  Figure  17  to  Figure  23);  increasing 
the  prf  decreases  the  viewing  area  (Figure  17  versus  Figure 
25);  decreasing  the  prf  increases  the  maximum  viewing  area 
(Figure  26  versus  Figure  17);  and  increasing  the  number  of 


Effects  of  Pulse  Width 


Figure  24  Effects  of  Pulse  Width  on  Blind  Zone 

prf  multiples  increases  the  maximum  range  (Figure  27  shows  a 
multiple  of  7  as  compared  to  Figure  28  which  shows  the  effects 
of  a  single  prf ) . 

Doppler  Maps .  The  relative  motions  of  the  transmitter, 
receiver,  and  target  dictate  the  doppler  shift  at  any  given 
point  in  X,  Y,  Z  space.  A  good  way  to  demonstrate  this  is 
to  view  a  plane  of  interest  and  see  the  doppler  contours  on 
that  plane.  For  example,  assume  the  ground  plane  is  the 
plane  of  interest.  The  doppler  shift  of  any  X,  Y  point 

.{Z  =  0}.  is  a  function  pf_the  velocity  of  _th^^  _ 

with  respect  to  the  transmitter  and  receiver  platforms. 
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Figure  27  Minimum  and  Maximum  Isoranges 
(Multiple  PRFj  NPRF  =  7) 


Pigure  28  “  Minimum  and  NaxiTnum  "Isoranges 
(single  PRFi  NPRF  =  1) 


The  doppler  map  developed  by  the  model  is  a  grid  of 
numbers  representing  the  doppler  shift  at  that  particular 
X,  Y  point.  Because  the  doppler  shift  is  well  behaved,  a 
contour  can  be  drawn  by  connecting  interpolated  points  on 
the  map.  Figure  29  shows  the  doppler  contours  for  a  receiver 
headed  to  the  left  at  600  KPH  and  a  transmitter  headed  to 
the  right  at  600  KPH.  The  relative  separation  is  1200  KPH. 
Figure  30  shows  the  contours  for  receiver  and  transmitter 
headed  left  at  600  KPH  and  remaining  a  constant  distance 
apart.  Figure  31  shows  the  same  pair  closing  with  equal 
velocity.  Figure  32  shows  the  transmitter  receiver  pair 
flying  orthogonal  (90  degrees  with  respect  to  each  other). 
Note  the  skewing  of  the  doppler  contours . 


Figure  29  Doppler  for  Transmitter"'and  Receiver  'Sepafating 
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Figure  32  Doppler  for  Transmitter  and  Receiver  Orthogonal 

The  transmitter  receiver  pair  will  not  always  fly  such 
regular  patterns,  nor  will  they  always  have  the  same  veloc¬ 
ity.  Figure  33  shows  the  effects  of  equal  velocities  but 
arbitrary  headings.  Figure  34  shows  the  effects  of  those 
same  headings  but  with  different  receiver  and  transmitter 
velocities . 

For  the  case  of  a  grovjnd  transmitter,  the  velocity  of 
the  transmitter  is  the  same  as  the  ground  plane.  The  trans¬ 
mitter  is  at  an  altitude  of  10  meteors.  The  doppler  shift  is 
now  only  a  function  of  the  receiver  velocity  and  heading. 
Figure  35  shows  the  doppler  map  for  a  receiver  separating  at 
600  KPH,  Figure  36  shows  the  receiver  closing  on  the  trans¬ 
mitter,  Note  the  change  in  sign  of  the  contours. 
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Figure  37  shows  the  contours  for  a  rei...  flying 
orthogonal  to  the  transmitter.  Figure  38  d  .  .  Jtrates  the 
skewing  of  the  doppier  contours  when  the  receiver  is  at  a 
heading  45  degrees  with  respect  to  the  receiver  transmitter 
baseline. 


Figure  37  Peceiver  Orthogonal  to  Ground  Transmitter 
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Figure  38  Receiver  45  Degrees  from  Ground  Transmitter 

Signal  to  Clutter  Calculations 

In  Chapter  III  th*'  Signal  to  Clutter  Ratio  (SCR)  was 
developed.  The  model  parameters  were  varied  to  demonstrate 
the  effects  of  changing  radar  cross  section  and  beam  sizes. 
Figure  39  shows  the  signal  to  clutter  ratio  for  the  generic 
waveform.  The  elevation  angle  is  displayed  in  5  degree 
increments.  SCR  is  shown  as  a  function  of  azimuth  and  eleva¬ 
tion.  Figure  40  demonstrates  the  change  in  the  SCR  when  a 
large  target  is  viewed.  Table  II  lists  values  for  the  SCR. 
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TABLE  II 


Signal  to  Clutter  Values  Generic  Waveform 


Azimuth  Angle 


-180 

-150 

-120 

-90 

-60 

-30 

0 

30 

60 

90 

120 

150 

180 

El 

-37. 

-35. 

-33. 

-30. 

-23. 

-42. 

-55. 

-42. 

-23. 

-30. 

-33. 

-35. 

-37. 

5 

-36. 

-34. 

-31. 

-26. 

-33. 

-45. 

-53. 

-45. 

-33. 

-26. 

-31. 

-34. 

-36. 

10 

-36. 

-33. 

-30. 

-24. 

-35. 

-44. 

-51. 

-44. 

-35. 

-24. 

-30. 

-33. 

-36. 

15 

-36. 

-33. 

-29. 

-23. 

-35. 

-44. 

-50. 

-44. 

-35. 

-23. 

-29. 

-33. 

-36. 

20 

-36  . 

-33. 

-29. 

-22. 

-35. 

-43. 

-50. 

-43. 

-35. 

-22. 

-29. 

-33. 

-36. 

25 

-36. 

-33. 

-28. 

-22. 

-35. 

-43. 

-49. 

-43. 

-35 . 

-22. 

-28. 

-33. 

-36. 

30 

-36. 

-32. 

-28. 

-21. 

-34. 

-42. 

-48. 

-42. 

-34. 

-21. 

-28. 

-32. 

-36. 

35 

The  use  of  a  ground  transmitter  changes  the  SCR  slightly. 
Figure  41  shows  the  SCR  for  a  small  target,  ground  transmit¬ 
ter  combination.  Figure  42  shows  the  SCR  for  a  large  tar¬ 
get.  Table  III  lists  the  values  for  tho  generic  waveform 
and  a  ground  transmitter. 


Ficurc  41 


Ground  Transmitter,  Small  Target 
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Figure  42  SCR;  Ground  Transmitter,  Target  RCS  =100  Sq.  M. 

TABLE  III 


SCR  Values  Generic  Waveform,  Ground  Transmitter 


Azimith 

Angle 

> 

180 

-150 

-120 

-90 

-60 

-30 

0 

30 

60 

90 

120 

150 

180 

El 

37. 

-35. 

-33. 

-30. 

-23. 

-42. 

-55. 

-42. 

-23. 

-30. 

-33. 

-35. 

-37. 

5 

36. 

-34. 

-31. 

-26  . 

-33. 

-45. 

-53. 

-45. 

-33. 

-26. 

-31. 

-34. 

-36. 

10 

36. 

-33. 

-30. 

-24. 

-35. 

-44. 

-51. 

-44. 

-35. 

-24. 

-30. 

-33. 

-36. 

15 

36. 

-33. 

-29. 

-23. 

-35. 

-44. 

-50. 

-44. 

-35. 

-23. 

-29. 

-33. 

-36. 

20 

36. 

-33. 

-29. 

-22. 

-35. 

-43. 

-50. 

-43. 

-35. 

-22. 

-29. 

-33. 

-36. 

25 

36. 

-33. 

-28. 

-22. 

-35. 

-43. 

-49. 

-43. 

-35. 

-22. 

-28. 

-33. 

-36. 

30 

36. 

-32. 

-28. 

-21. 

-34. 

-42. 

-48. 

-42. 

-34. 

-21. 

-28. 

-32. 

-36. 

35 

Doppler  Spread  of  the  Clutter  Cell 

The  doppler  spread  of  the  clutter  cell  is  a  function  of 
the  beam  shape  and  the  _platform  velocities  .  The  clutter,  cell 
doppler  spread  for  a  separating  receiver  transmitter  pair 
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The  doppler  spread  varies  with  the  flight  path  orienta¬ 
tion.  If  the  transmitter  and  receiver  are  flying  orthogonal 


the  doppler  spread  skews  as  shown  in  Figure  46,  Figure  47 
shows  the  doppler  spread  for  a  transmitter  receiver  pair 
flying  arbitrary  headings  of  53  degrees  and  315  degrees 
respectively.  Figure  48  shows  the  doppler  spread  for  trans¬ 
mitter  and  receiver  having  the  arbitrary  headings  and  dif¬ 
ferent  velocities. 

DOPPLER  SPREAD  OF  BEAil  (AZIMUTH  AHGLE  US  ELEUATIOh  ANGLE) 


Figure  46  Doppler  Spread,  Orthogonal 


The  doppler  spread  for  a  ground  transmitter  to  airborne 
receiver  is  dependent  upon  the  beam  shapes  of  both  platforms 
but  only  the  velocity  and  heading  of  the  receiver.  There¬ 
fore,  the  doppler  spread  of  the  air  to  ground  system  is 
better  behaved.  Figure  49  shows  the  doppler  spread  for  a 
closing  receiver  while  Figure  50  shows  the  doppler  spread  for 
a  separating  pair.  Note  that  a  closing  and  separating  sys¬ 
tem  have  the  same  spread.  It  should  be  realized  that  the 
doppler  spread  is  the  same  but  the  doppler  shift  of  the  clut¬ 
ter  cell  is  not  at  the  same  place  in  the  frequency  spectrum. 
If  the  receiver  flys  orthogonal  to  the  transmitter,  the  dop¬ 
pler  spread  of  the  beam  shifts  as  shown  in  Figure  51. 


Figure  49  Doppler  Spread,  Air  to  Ground,  Closing 


Tlie  aii'  to  ground  ^i/steir.  doppler  spread  shifts  with 
respect  to  the  heading  of  the  receiver.  Figure  52  shows  the 
shift  for  a  heading  of  3l5  degrees  with  respect  to  the  trans 
mitter  receiver  baseline  'hile  Figure  53  shows  the  doppler 
shift  for  a  heading  of  243  degrees. 

I'OPPLES  '•PhFhL  op  bEPil  vrtZIfiUTH  ft'rlLE  OS  ELEOATION  ANGLE) 


TABLE  IV 


Doppler  Shift  Versus  Velocity 


Doppler 
Frequency 
Shift  (Hz) 

Kilometers 

Per  Hour 

Miles 

Per  Hour 

Nautical  Miles 
Per  Hour 

2000 

720 

446 

389 

1500 

540 

335 

292 

1000 

360 

223 

194 

500 

180 

112 

97 

269 

97 

60 

52 

Sample  Scenario 

Suppose  that  the  doppler  spread  and  the  bistatic  detec 
tion  regions  for  a  possible  surveillance  scenario  are  of 
interest.  Assume  two  similar  radar  platform  headed  in 
opposite  directions,  passing  side  by  side  separated  by  175 
kilometers  when  orthogonal  as  shown  in  Figure  54. 


Figure  54  Possible  Air  to  Air  Scenario 


Figure  55  is  a  composite  of  three  bistatic  detection 
region  graphs  to  show  the  area  of  coverage  from  the  time  the 
base  line  is  200  Kilometers  until  the  platforms  are  at  right 
angles  separated  by  175  Kilometers.  As  the  platforms  sep¬ 
arate,  a  mirror  image  (with  respect  to  the  graph  base  line) 
of  the  detection  region  is  traced  out.  The  time  between 
curve  one  and  curve  two  is  261  seconds.  The  time  betvreen 
curve  two  and  curve  three  is  30  seconds.  Note  changes  in 
both  areas  which  are  visible  and  blind  to  the  receiver. 


Jusl  as  the  detection  regions  change,  the  doppler 
spread  of  the  clutter  ceii  changes.  Figure  56  shows  the 
clutter  cell  doppler  spread  -for  initial  acquisition  at  200 


kilometers.  Figure  57  shows  the  change  in  ;.he  doppler  spread 
as  the  transmitter  receiver  pair  approach  their  minimum 
separation,  while  Figure  58  shows  the  doppler  spread  during 
the  moment  when  the  pair  are  orthogonal.  Figure  59  shows  the 
doppler  spread  as  the  transmitter  and  receiver  platforms 
separate.  Note  in  particular  the  doppler  shift  for  the  var¬ 
ious  angles  of  beam  elevation  and  the  shift  in  the  magni¬ 
tude  of  the  beam  spread. 

\ 

DOPPLER  SPREAD  OP  BEAM  CAZlMUtH  ANGLE  US  ElEUATlON  ANCLE) 
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Figure  56  Doppl-er  Spread,  Closing 
(Separation .=  200  KM) 


DOPPLEk  spread  of  beam  (AElflUTH  ANGLE  OS  ELEOATION  ANGLE) 


The  sample  scenario  shows  the  shift  in  the  detection 
regions  and  the  shift  in  the  doppler  spread  of  the  clutter 
cell.  A  passive  receiver  must  be  capable  of  dynamically 
adjusting  parameters  in  order  to  maximize  the  possible  infor¬ 
mation  available  at  its  antenna. 


V.  Summary.  Conclusions  and  Recommendations 


Summary 

The  bistatic  radar  problem  is  composed  of  many  varia¬ 
bles.  The  model  developed  demonstrates  the  need  for  under¬ 
standing  the  effects  of  these  variables  on  the  detection 
regions.  The  model  demonstrates  the  ability  to  predict  the 
signal  to  clutter  ratios  and  doppler  effect  for  particular 
beam  shapes  and  system  configurations.  The  system  and  its 
viewing  area  is  more  complex  and  limited  when  compared  to  a 
monostatic  radar  system.  The  bistatic  system  is  faced  with 
a  blind  zone  which  is  controlled  by  the  waveform  parameters. 
Chapter  IV  examples  demonstrated  the  effects  of  varying  the 
parameters  such  as  pulse  duration  and  pulse  repetition  fre¬ 
quency.  Varying  other  system  parameters  such  as  the  antenna 
gain  and  beam  shape  have  similar  effects  on  the  Ovals  of 
Cassini  . 

The  figures  of  merit  for  clutter  and  doppler  spread  of 
the  beam  are  predictable,  and  therefore  usable  in  real  time 
processing.  The  doppler  spread  and  SCR  functions  are  well 
behaved  and  will  allow  a  surveillance  receiver  to  continually 
update  the  waveform  pr<x:essing  portion  of  the  receiver,  Cwi- 
tinuous  updates  would  allow  for  improved  performance  and 
sub-clutter  visibility  by  having  a  dynamically  varying  clut¬ 
ter  notch  filtei. 

The  means  for  exploiting  a  bistatic  radar  system  are 
apparent.  If. is  system  is  most  vulncralle  alojio  the 
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transmitter  receiver  baseline.  The  blind  zone  is  limited 
by  waveform  and  doppler  spread.  Additionally,  the  system 
needs  to  be  able  to  have  rapid  synchronization  and  acquisi¬ 
tion  of  the  transmitter  waveform,  changing  the  waveform 
will  disrupt  the  receiver  surveillance  until  the  system  can 
be  resynchronized.  Any  parameter  which  maximizes  the  viewing 
area  of  the  receiver  platform  is  a  prime  candidate  for  change 
if  the  transmitter  wishes  to  deny  or  decrease  the  viewing 
area  of  the  receiver.  Specifically,  increasing  the  pulse 
width  or  the  prf  and  decreasing  the  prf  multiple  will  cause 
a  drastic  change  in  the  receivers  viewing  area.  In  fact, 
these  changes  could  effectively  blind  the  receiver. 

Conclusions 

Choice  of  a  good  host  receiver  pair  waveform  is  essen¬ 
tial  in  maximizing  the  bistatic  viewing  area.  For  example, 
it  was  shown  in  Chapter  IV  that  a  low  prf,  multiple  prf, 
short  pulse  duration  waveform  provides  much  greater  viewing 
area.  Examination  of  the  equations  of  interest  lead  to  the 
conclusion  that  to  maximize  the  viewing  area,  the  system 
must  strive  to  maximize?  power,  gain  of  the  antennas  and 
the  bistatic  radar  cross  section  of  the  target .  The  system 
should  also  minimize?  the  receiver  bandwidth,  the  transmitted 
wave  length  {maximize  transmitter  frequency),  system  losses, 
and  receiver  front  end  temperature. 

The  ciioice  of  maximizing  or  nu?umizing  a  particular 
parameter  may  not  be  within  the  capability  of  the  r^eiver 
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if  the  transmitter  is  non-cooperative.  Likewise,  the 
receiver  can  do  nothing  to  change  the  radar  cross  section 
of  a  target.  The  practical  step  is  to  maximize  those  para¬ 
meters  which  are  within  the  control  of  the  receiver.  If  the 
transmitter  and  receiver  are  cooperative,  the  transmitter 
must  generate  a  waveform  which  is  beneficial  to  the  receiver. 
If  the  transmitter  is  non-cooperative,  the  receiver  should 
look  to  acquire  only  those  waveforms  which  provide  adequate 
coverage.  In  the  event  of  an  inadequate  transmitter 
receiver  pair  waveform,  the  receiver  would  have  to  generate 
its  own  waveform  and  become  a  monostatic  platform,  or  be 
prepared  to  cancel  its  surveillance  activities. 


Recommendations 

1.  The  Bistatic  Model  should  be  extended  to  include! 

-  The  Curved  Earth  Model 

-  Varying  Terrain  Backgi'ounds 

-  Sky  Clutter 

-  Varying  Beam  Shapes 

-  Side  Lobe  Detection  Regions 

-  Side  Lobe  Clutter 

-  Doppler  Range  Gate  Limitations  on  clutter 

-  Varying  Clutter  for  Different  frequencies 

2.  A  bistatic  radar  cross  section  model  should  be 
included  in  the  model.  The  cross  section  model  should  be  a 
fund  ion  of  look  elevation  angle  and  azimuth  angle.  The 
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Appendix  At  Fortran  Program  Listing 


C  tttttttttttttttttttttttttttttttttttttttttttttUtttttttttttttttt 

c 

C  BISTATIC  RADAR  PROGRAM 

C 

c 

1  DIMENSION  TABLE <300, 250) ,SPD( 10,37) , BASE (302) , ALT ;252>, 

C  BASEl <38) ,ALT2<38) ,11 <30), 12(63) , 13(62) , lyy (37) ,YY(37> 

C 

C  THE  DATA  statements  set  up  the  plot  10  labels 
C 

data  11/9,9,9,9,66,73,83,84,65,84,73,67,9,68,69,84,69, 

C67 ,84,73,79*78,9,82,69,71,73,79,78, 83/ 
data  12/9,9,9,9,9,83,73,71,78,65,76,9,84,79,9,67,76,85, 

C84, 84, 69, 82, 9, 82, 65, 84, 73, 79, 9, 40, 65, 90, 73, 77, 85, 84, 72, 
c9, 65, 78, 7 1,76, 69, 9, 86, 83, 9, 69, 76, 69, 86, 65, 84, 73, 79, 78, 
c9, 65, 78, 71, 76, 69, 41/ 

data  13/9,9,9,9,9,68,79,80,80,76,69,82,9,83,80,82,69,65, 

C68, 9, 79, 70, 9, 66, 69, 65, 77, 9, 40, 65, 90, 73, 77, 85, 84, 72, 
c9, 65, 78, 71, 76, 69, 9, 86, 83, 9, 69, 76, 69, 86, 65, 84, 73, 79, 78, 
c9, 65, 78, 71, 76, 69, 41  / 

2  PEAL  IS0R,IS0MIN, ISOMAX, LA, LN, LOSS 

3  INTEGER  NPRF,PRF 
c 

C  Initialize  the  free  variables  to  default  values 

c 

•'  P=O.01 

61=^  ?0 

6  GR=*i‘C» 

7  SNR=6.5 

8  RCS*1.75 

9  BW«=1 

10  L0SS«6.5 

15  TEMP*300 

25  XWAZ>1.5 

30  RWAZ-1.5 

31  RWEL-6 

32  XVCL-6 

33  F0«3 

34  Z=0 

35  VX=600 

36  XAZ=0 

37  XEL=^0 

38  VT«600 

39  TAZ=0  ‘  ’  ■ 

40  TEL=0 

41  TPULSE=5 


el 

-  't| 

f,| 


PRF=2SOOO 

NPRF«3 

XALT-9 

RALT»9 

SEP=200 

C=3E5 


dummy4  is  a  parameter  that  the  user  sets  to  alloM  the  progras 

to  set  itself  up  for  graphics  or  Crt  mode 

dummyA  »  1  seids  plotlO  codes  to  the  tektronix  terminal 

PRINT*, 'ENTER  1  IF  YOU  ARE  ON  A  GRAPHIC  TERMINAL,  ANY  OTHER 
C  NUMBER  FOR  A  CRT’ 

READ* , DUMMY* 

call  welcoffl<dummy4) 

if  (dummy4. eq. 1 )  then 

t******************************************************** 

set  up  for  plotlO  graphs 

t*******!************************************************* 

ALT (252) *125 
ALT(2)=-125 
BASE(l)=301 
ALT(1)=251 
ALT2 ( 1 ) *37 
BASEl (l)=37 
DO  54  1*2,38 
BASEl (I)=-200+I»10 
CONTINUE 
DO  55  1*2.302 
BASE(I)*l-52 
CONTINUE 
DO  57  1*3,251 
ALT{I>*0 
CONTINUE 
end  if 

****** ************************************** t******************* 
The  user  selects  the  option  for  the  program  to  execute 

•ttSt*«**«*»***»**«*C****«****S**S*t*«**«t*tt*tttt«tt**»tttt»ttt 


PRINT*, 'DO  YOU  WANTt  1 

print*, •  2 

PRINT*,'  3 

print*, '  4 

READ* , DUMMY3 
IF  (DUMMY3.BT,4>  THEM 
CALL  ERROR 
60  TO  64 

ELSE  IF  (DUMMY3.lt. I )  THEM 
CALL  ERROR 
GO  TO  64 
END  If 


1  ISORANGCS* 

2  tOPPLim  HAP' 

3  SIGNAL  TO  CLUTTER  CALCULATION' 

4  DOPPLER  SPREAD  OF  BEAM* 


200  i-f  <dummy4.  eq.  I )  then 

call  wipe 
end  if 

C 

C  The  program  prints  out  the  current  parameter . values  befor  any 
C  execution.  This  allows  the  operator  to  change  any  value  or 
C  set  of  values. 

C 

201  PR 1 NT ’CURRENT  VALUES  ARE  DEFINED  ASi’ 

202  PRINT!,’  TARGET  TRANSMITTER 

C  RECEIVER’ 

208  PRINT  1110, VT,VX, TEMP 

210  PRINT  1112.TAZ,XAZ,SNR 

214  PRINT  1114,TEL.XEL 

216  PRINT  1116,RCS,P,L0SS 

220  PRINT  1120,Z,XALT.SEP 

230  PRINT  1130,GT,6R 

240  PRINT  1140,F0,BW 

245  PRINT  1150, XWAZ,RWAZ 

247  PRINT  1118,XWEL,RWtL 

249  PRINT  11 19,PRF,RALT 

PRINT  llZl.TPULSE 
PRINT  1122,NPRF 

250  PRINT!, ’ENTER  OPTION  NUMBER  tl-26D  TO  CHANGE  VALUE’ 

260  PRINT!, ’OR  27  TO  CONTINUE* 

270  PRINT#, 'OR  28  TO  ABORT’ 

C  !!!!!!!!!!!!!!!!!#!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 

C 

C  DuiTuny  directs  the  program  to  an  execute  mode,  a  change  parameter 
C  mode  or  a  quit  mode 
C 

C  !!!!!!!!!!!!!#!!!!!!!!!!!!!!!!!!!!!!*!!!!!!!!!!!!!!!!«!! 

300  READ!. DUMMY 

310  IF  (DUMMY  .EG,  28)  THEN 

320  60  TO  7777 

330  ELSE  IF  (DUMMY  .EO.  27)  THEN 

340  60  TO  595 

350  ELSE  IF  (DUMMY  .6T.  28)  THEN 

355  CALL  ERROR 

360  GO  TO  250 

370  ELSE  IF  (DUMMY  .LT.  I)  THEN 

380  CALL  ERROR 

390  60  TC  250 

400  END  IF 

c 

c  subroutine  call  to  change  parameters 

c 

50C»  CALL  CHANGE  (DUMMY, VT.TAZ,  TEL,  VX,  XAZ,  XEL,TPULSE,PRF,NPRF, 
cFO ,  SEP ,  Z ,  X ALT ,  F ,  GT ,  (iR  .  SNR ,  RCS ,  BW ,  L OSS ,  TEMP ,  RWEL  .  RWAZ  , 

CXWEL. XWA:,RmL'  > 


c  rvcycl*  to  print  out  oi  cuurcnt  value  table 
c 

505  GO  TO  200 

C  TEST  DUMMY3  TO  FIND  PROGRAM  TO  RUN 
595  IF  (DUMMY4.Eq. 1)  THEN 
CALL  INITT<1200) 

CALL  BINITT 
end  if 

C 

C  vector  to  program  selection  based  on  dumayS 
C 

597  go  to  (600,640,744,744)  DuiiUBy3 

C 

C  subroutine  call  to  calculation  power  equation  limitation 

C  value,  R4 

C 

600  CALL  RANGE (P,6T,6R, SNR, RCS,FO,LOSS,BW, TEMP, R4) 
c 

c  subroutine  call  to  determine  the  time  domain  limitations 

C  of  the  minimum  and  maximum  isoranges 

c 

610  CALL  MINMAX <TPULSE,C,SEP,NPRF,PRF, ISOMIN, ISOHAX, I.XALT.RALT) 

c 

c  if  not  on  a  40XX,  print  values  for  isooin  and  isomax  on  screen 
c 

611  if (dummy^. ne. 1 )  then 

print  1160,isomax 
print  1161,isoffiin 
go  to  990 
end  i  f 
c 

c  if  on  a  40XX,  draw  the  power  equation  limitation  regions 
c  (  the  Owals  of  Cassini) 

c 

612  call  check (base, alt) 
call  dsplay (base,alt) 
call  hi abel (30,11) 

615  CALL  POWER (Z,X ALT, SEP, R4, BASE, RALT) 
c 

c  subroutine  call  to  draw  the  isorange  ellipsiods 

c 

620  CALL  ISORdSOMIN,  ISOMAX,  1, SEP,  Z,  XALT. BASE, RALT) 
c 

c  a  check  to  see  if  all  the  intermediate  ellipsiods  are  drawn 

c  if  not,  the  program  will  recycle  through  the  isor  routine 

c  until  the  minimun,  maximum  and  intermediate  ellipses  have 

c  been  dravm 

c 

617 


IF  d.iiu.Oi  THEN 


630 

635 

636 
640 


c 

c 

c 

c 

645 

c 

c 

c 

c 

c 

646 

647 
C 

C 

c 

C 

c 

744 


745 

750 

755 

C 

c 

c 

c 

775 


C 

C 

c 

c 

780 

c 

C 

c 

785 

790 


GO  TO  620 
END  IF 
go  to  990 

IF  (DUI1MY4.EQ.  DTHEN 
call  Mi pc 

END  IF 

subroutine  call  to  calculate  the  dc^pler  values  for  the  region 
defined  in  the  model. 

CALL  DOPMAP ( VX, XEL. XAZ, VT, TEL, TAZ, FO, C, SEP, XALT, 

CZ, TABLE, RALT) 

a  subroutine  call  ♦‘o  print  the  doppler  values  for  the  plane  of 
interest,  if  a  graphical  plot  package  is  added  at  a  later  date, 
it  should  be  added  here. 

CALL  DOPPLT (TABLE) 

60  TO  990 

THE  SCR  ROUTINE.  This  routine  calculates  the  SCR  for  seven  look 
down  angles  (5,10,15,20,25,30,35).  The  scr  is  a  function  of  the 
Azimuthal  scan  angle  as  well. 

RCSDB= I Ot LOG 1 0 (RCS ) 

IF(Z.GE.RALT)  THEN 

PRINTt,'  MODEL  HAS  NO  CAPABILITY  TO  FIGURE  "SKY“  CLUTTER* 

60  TO  990 
END  IF 
DO  900  1-1,7 
RSEL-It5 

TABLE (1 ,38 ieRSEL 
DO  800  J*l,37 

RSAZ»(J-*1)  •  10-180 
Iyy(J)«INT(RSAZ> 

The  subroutine  call  to  calculate  the  doppler  spread  of  the 
cocaon  beam  area. 

CALL  DOPSD(S€P,2.XALT,RSAZ.RSEL,RllA2.RICL,XWAZ,XlieL, 

C  VX,VT,FO,C,I)OPSR,XS£L,XRAZ,PMXEgU,XEL,XAZ,TA2,T€L,PHI,R;^T> 
SP0(X,J)-B0PSR 

subroutine  call  to  calculate  the  per  unit  density  of  the 
common  beam  area 

CALL  DENT (PHI, PHIEOU, DEN) 

subroutine  call  to  calculate  the  common  beam  area  -  •  - - 

CALL  ARr«  (SEP,  Z,P^L,RrA7,PW,  XWAZ.RWAr  ,C,  ACDP.RALT,  XALT) 
table ( 1 , J ) -RCSDB-DEN-ACDB 


c  ^ 


800  CONTINUE 

900  CONTINUE 
c 

C  This  section  sets  up  the  plotlO  call  for  drawing  the  SCR  graph 
c 

901  IF  (DUMMYA.Eq.l)  THEN 

if  {dummy3.eq.3>  then 
DO  902  I«=l,37 

ALT2(I+l)»table(l, I) 

902  CONTINUE 

CALL  CHECK ( BASE 1,ALT2) 

CALL  DSPLA Y  <  EASE 1 , ALT2 ) 
call  hlabel (63,i2) 
do  904  i=2,7 

do  903  i«l,37 

alt2<I+l)=TABLE<J,i) 

903  continue 

call  cplot (basel ,alt2) 
n=int ( j/2) 
call  llne(n) 

904  continue 

go  to  990 
end  if 

DO  905  1=1,37 
ALT2(I  +  l)=spd(l,  1) 

905  CONTINUE 

call  hlabel (62,  i3) 

CALL  LINE (34) 

CALL  CHECK < BASE 1,ALT2) 

CALL  DSPLAY ( BASE 1 , ALT2 ) 

DO  907  1=2,7 

DO  906  J=l,37 

ALT2(J+l>=spd<I, J) 

906  CONTINUE 

CALL  CPL0T(BASE1,ALT2) 

907  CONTINUe 

po  to  990 
end  if 
C 

Z  Tf  not  on  a  40xXy  the  program  will  Print  a  table  of  SCR 
c  values  and  dopplcr  values  on  the  CRT  screen. 

C 

910  PRINT*,!'  SIGNAL  TO  CLUTTER  RATIO  AND  DOPPLER 

C  SPREAD') 

PRINT  1162, <Iyy<I) ,1=1,37,3) 

915  DO  955  1=1,7 

920  PRINT  1000, (TABLE(I,J),J=1, 37,3) ,TABLE(I, 38) 

925  PRINT  lOCl, (SPD(I,J) ,3=1,37,3) 

^55  ■  -CONI  INUE . - . -  - - - - - - 

990  if  (dummy4.eq. 1)  then 
call  tinput<y) 
call  finitt (0,700) 


end  if 

PRINTt,*DO  YOU  WISH  TO  CHANGE  PARAMETERS  Ai4D  TRY  AGAIN 
C  <ENTER  1)’ 


995 

READ*, DUMMY 

996 

IF  (DUMMY. EO.l)  THEN 

997 

GO  TO  60 

998 

END  IF 

1000 

FORMAT (»  SCR  dB  »14F5.0) 

1001 

FORMAT <*  SPREAD  »,13F5.0) 

1110 

FORMAT <»  1 

VELOCITY  * 

’,F7.2, ’ 

6 

VELOCITY 

CF7.2,'  18 

TEMP  »= 

’ ,F6.2) 

1112 

FORMAT!’  2 

A2  ANGLE  * 

’,F7.2,’ 

7 

AZ  ANGLE 

m 

CF7.2,’  19 

8NR  « 

’ ,F6.2) 

1114 

FORMAT!’  3 

EL  ANGLE  = 

’,F7.2,’ 

8 

EL  ANGLE 

» 

CF7.2> 

1116 

FORMAT!’  4 

RCS  = 

’,F7.2,’ 

9 

POWER 

m 

CF7.4,’  20 

LOSS  = 

’,F7.2> 

1120 

FORMAT!’  5 

2  PLANE  « 

’ .F7.2,’ 

lO 

ALTITUDE 

B 

CF7.2,’  21 

SEP  = 

’,F7.2) 

1130 

FORMAT ( ’ 

11 

GAIN 

B 

CF7.2,’  2 

GAIN  = 

’,F7.2) 

1140 

FORMAT 

12 

CARRIER 

t= 

CF7.2,’  23 

BW  « 

’,F7.2) 

1150 

FORMAT ( ’ 

13 

BEAM  WD 

B 

CF7.2,’  24 

BEAM  ND  = 

’,F7.2) 

1118 

FORMAT !  ’ 

14 

BEAM  HT 

CF7.2,’  25 

BEAM  HT  = 

’,F7.2) 

1119 

FORMAT ! ’ 

15 

PRF 

ss 

CI7,’  26 

ALT  <=  ’ , 

F7.2) 

1121 

FORMAT ! ’ 

16 

PULSE  WD 

s 

CF7.2) 

1122 

FORMAT !  ’ 

17 

«  PRF 

B 

CI7) 

1160 

format!’  Maximum  Isorange  «  ’,f6. 

2) 

1161 

format!’  Minimum  Isorange  *  *,f6. 

2) 

1162 

f  ormat ! ’  Az 

’,1315,* 

El  ’) 

7777 

PRINT*, 'THANK  YOU  AND  COME  AGAIN’ 

STOP 

END 


c 

c  e  short  routine  to  send  «n  error  eesssQe  to  the  operstsr 
c  if  an  out  of  bounds  parameter  is  specified. 


c 

1 

SUBROUTINE  ERROR 

5 

PRINT*, ’ERROR  IN  DATA  ENTRY, 

TRY  AGAIN’ 

10 

RETURN 

15 

END 

c 

£  -a 

short  -welcome 

c 

1  SUBROUTINE  WELCOM (dummyA) 

PRINT*,'  ' 


_  ejc  _ 


PRINT* » 

PRINT*,*  Welcome  to  the  World  of  Dietetic  Redar* 

PRINT*,*  * 

PRINT*,*  » 

PRINT*, ’The  purpose  of  this  program  is  to  assist  the  study* 
PRINT*, *of  Bistatic  Radar.  This  program  allows  for  varying* 
PRINT*, ’parameters  to  develop  isodoppler  and  isorange  owps,* 
print*, *as  well  as  Signal  to  Clutter  Ratio  values  and  the  ,* 
print*, ’doppler  frequency  spread  of  the  clutter  cell.* 

PRINT*, *The  maps  can  be  displayed  on  a  Tektronix  40XX  graphics* 
PRINT*, ’terminal  and/or  generated  in  hardcopy.  The  program* 
PRINT*, 'will  ask  you  to  input  data.  Enter  the  data  required:* 
PRINT*, ’All  distance  moasures  are  in  Kilometers.  All  angles’ 
PRINT*, ’in  Degrees.  All  Velocities  in  Kilometers  per  Second.* 
PRINT*,’  ' 

PRINT*, ’The  Bi static  Radar  system  modeled  by  this  system  always 
PRINT*, ' assumes  the  receiver  and  transmitter  form  the  baseline. 
PRINT*, ’All  angles  are  measured  with  reference  to  the  base  line 
PRINT*,’  ’ 

PRINT*,’  /» 

PRINT*,’  /  Angle  Theta  (positive)* 

PRINT*,’  /* 

PRINT*,  ’  **»*«***»-i***»**|******»******»»**«******«*** 

PRINT*, ’ Receiver  Transmitter* 

PRINT*,’  ’ 

PRINT*, ’STRIKE  ANY  NUMBER  AND  CARRIAGE  RETURN  KEY  TO  CONTINUE* 
READ  *, DUMMY 
if  (dummy4, eq. 1 )  then 
call  wipe 
end  i  f 
PRINT*,’  ’ 

PRINT*,’  ’ 

PRINT*, ’Remember: ’ 

PRINT*,*  * 

PRINT*, ’The  relationship  of  kilometers  to  nautical  miles  isi' 
PRINT*, ’I  Nautical  Mile  =  1.852  Kilometers* 

PRINT*, 'The  relationship  of  kilometers  to  statute  miles  isi* 
PRINT*,  *1  Statute  Mile  «  1.613  Kiloaieters* 

PRINT*, ’The  relationship  of  meters  per  second  to  nautical  miles 
PRINT*, *por  hour  isi  1  NM/Mr«  0.514  M/SEC* 

PRINT*, ’The  relationship  of  meters  per  second  to  statute  milms* 
PRINT*, ’per  hour  is:  1  M/Hr  »  0.448  H/GEC* 

PRINT*,’  * 

PRINT*,'  ’ 

RETURN 

END 

THIS  SUBROUTINE  CALCULATES  THE  DENSITY  OF  THE  CLUTTER  CELL 


SUBROUTINE  DENT (PHI .PHIEOU, DEN) 
X«PHI 

IF  (X.GT.VO)  THEN 


SIS  DEN«SQRT(<X-90)*3.2111)-32 

40  ELSE 

4b  DEN*SQRT ( (90-X) »10) -32 

EO  END  IF 

S5  IF(PHIEQU.6T.30)  THEN 

5A  PRINT*, ’MODEL  ILL  DEFINED  AT  THIS  POINT' 

CORR-0 

57  ELSE  IF <PHIEQU.lt. 1.5)  THEN 

58  C0RR=<PHI/180)*(-8) 

59  ELSE 

CORR= ( PH I / 1 80  >  * ( 9 . 2*L0G 1 0 ( PH I EQU) -9.2) 

END  IF 

DEN  >  DEN  +  CORR 

60  RETURN 

65  END 

C 

C  THIS  SUBRCUTINE  CALCULATES  THE  Area  of  the  clutter  cell 
c  which  is  Jependent  upon  the  bean  geometry  and  the 
c  I pok  down  and  az i muthal  angl  es 

c 

1  SUBROUTINE  AREA (SEP, Z, RSEL, RSAZ , BW, XWAZ , RWAZ , C, ACDB, RALT, XALT) 

5  R2’  <RalT-Z)/(SIND  ^RCEL>) 

10  Y=R2»C0SD(RSEI  )*SIND(RSA7) 

15  X=R2*nOlO (RSEL) »COSD (RSAZ) 

20  R1=SQRT( (SEP-X)»12+V**2+(XALT-Z)»*2) 

22  rwaz  l^^rwaztO.OlTAS'-’iS 

23  1  ex  wt^z  *0.0 174533 

25  R3=SQRT((SEP-X)**2+<XAL1-RALT)**2) 

30  SEC2E:ie4»Rl*R2/(2tRi*R2+R2»*2+Rlt*2-R3**2) 

35  TEST  e  R2*RWAZt/ni 

AO  IF(TEST.GE.XWAZl)  THEN 

45  AC=Rl*XWAZl»C*SeC2Bl/<2*BW) 

50  ELSE 

55  AC=R2*RWAZ:tC*SEC2BI/<2*BW: 

60  END  TF 

100  ACODelO*LOG10(AC) 

105  RETURN 

no  END 

C  Conputc  The  Doppler  Spread,  and  the  Traneaitter 
C  Elevation  and  Azimuth  Angles 

1  SUBRCUTINE  DOPSD ( SEP, Z, XA.T, RSAZ, RSKiLtR^iAZyRNEL »%ifAZ, XNEO., 

CVX,  VT,FO,C,DOPSR,  XSE1.,XRAZ,PH1EQU,  XEL,  :<AZ,TAZ,t€L,rHt,RALn 

2  Dll^NSION  GATE<e,4) 

5  REAL  LA,LN 

25  RWEL2«RWEL*0.5 

30  RWAZ2eRWA2*0.5 

35  XWAZ2eXWA2*0.S 

40  XWEl.2=XWEL$0.5 

45  TEMP=(RALT-Z)/fAND<R^L') 

50  X»TeHP*COSD(RSAZ> 

55  V«T£MPiSIND<RSAZ) 

60  TEMF2eSQRT(Ytt2*(REP-X) *»2» 


101 


65 
C 

66 


67 

70 

C 

72 

75 

80 

85 

90 

95 

100 

105 

110 

115 

120 

125 

130 

135 

140 

145 

150 

155 

160 

165 

C 

c 

c 

170 

17S 

leo 

les 

190 

195 

200 

205 

210 

216 

217 

220 

2X1 

222 

225 

226 
-227 

230 

C 

235 


XSA2*ATAND ( Y/ (SEP-X> ) 

THE  AZIMUTHAL  ANGLE  18  NOW  DETERMINED 
IF  (RSAZ.6T.180>  THEN 
RSAZ=360-RSA2 
END  if 

PHI «ABS ( XSAZ  >  ^ABS (RSAZ ) 

XSEL»ATAND  <  <  XALT-Z ) /TEMP2) 

THE  EQUIVALENT  ANGLE  CAN  NOW  BE  DETERMINED 
PHIEQU- <PSEL+7SEL) /2 
6ATE  (1,1)  «^TEMP*COSD  (RSA2-RWAZ2) 

GATE  (2,1)  «TEMP*COSD  (RSAZ-fRWA22) 

GATE (3,1)* (RALT-Z ) »COSD (RSAZ ) »TAND (RSEL+RWEL2) 

GATE (4,1)* (RALT-Z ) tCOSD (RSAZ ) tTAN0(RSEL-RWEL2) 

GATE (1,2) *TEMP$S1ND (RSAZ-RWAZ2> 

GATE  (2, 2)  eTEJ1P»SIND  (RSA2+RWA22) 

GATE  (3,2)*  (RALT-Z  )  «SIND  (RSAZ  >  ITAND  (RSEL-I-RWEL2) 

GATE (4,2)* (RALT-Z ) tSIND (RSAZ ) ITAND (RSEL-RWEL2) 

GATE (5,1) *SEP-TEMP2*C0SD (XSA2-XWAZ2) 

GATE  <6, 1 } *SEP-TEMP2»C0SD (X8A2+XWAZ2) 

GATE (7.1 ) *SEP- <  XALT-2 ) $TAND ( XSEL-XWEL2) tCOSD ( XSAZ) 

GATE  (8,1)  *»SEP-  ( XALT-2 )  ITAND  ( XSEL4-XWEL2)  ICOSD ( XSAZ) 

GATE (5,2) *TEMP2*SIND (XSAZ~XWA22) 

GATE  (6,2)  *T£riP2ISIND  (XSA2+XWAZ2) 

GATE  <7 , 2)  =  ( XALT-Z )  ITAND  ( X13EL-XWEL2)  ISIND  (XWAZ2) 

GATE (8, 2) * ( XALT-Z ) ITAND ( XSEL+XWEL2) ISIND ( XWAZ2) 

do;  160  1^1,8 

GATE  ( 1 , 3)  =  ( X-6ATE  ( 1 , 1  > )  1 12-^  ( Y-6ATE  ( 1 , 2) )  112 
CONTINUE 

NOW  THE  PROGRAM  LOOKS  FOR  THE  SMALLEST  FOUR  DISTANCES  AND 
PLACES  THE  COORDINATES  IN  THE  TOP  FOUR  RUNGS  OF  GATE 
DO  230  N=l,7 

L=N 
JJ-N+l 

DO  210  I“JJ,8 

IF(6ATE(l,3>.LT.6ATE(I,3))  THEN 
GO  TO  230 
END  IF 
L*I 

CONTINUE 

T*GATE(L,3) 

T2*GATE(L,2) 

T3»6ATE(L,1) 

GATE  (L ,  3 )  *GATE  (N,  3 ) 

GATE  (L,  2) -GATE  (N,  2) 

GATE(L,1)«6ATE(N,  1) 

GATE(N,3)e=T 

GATE(N,2)-T2 

- 6ATE(N,1W3 - — . .  —  - 

CONTINUE 

NOW  FIND  THE  DOPPLER  OF  THE  FOUR  POINTS 
VX2  e  (VXISIND(XEL> ) /3600 


-  KC  - 


240 

245 

250 

255 

260 

265 

270 

275 

276 

277 

278 

279 

280 
281 
282 
290 
295 
300 
305 
310 
320 
330 
c 


c 

c 

1 

5 

10 


30 

35 

40 

45 

50 

55 

70 

100 

105 

110 

115 

c 

c 

c 

,  -.1 .  .. 


VXY  »  (VXt(C0SD(XEL))«8XN0(XAZ))/3600 

VXX  «  <9Xt(C08D(XEL)}«C08D(XAZ}>/3600 

VTZ  «  (VTt8IND(TEL)>/3600 

VTY  *=  (VT«(C08D(TEL>)«SIND(TAZ})/3600 

VTX  «  (VT*<COSD(TEL))*COSD(TAZ)>/3600 

LA«C/(F0*1E9) 

DO  295  1-1,4 
X«6ATE<I,1) 

Y»GATE<I,2) 

AP  =  SQRKXtX  ♦  Y*Y  ♦  <Z-RALT)  t  (Z-RALT) ) 

BP  -  SORT  ( (SEP-X)»(SEP-X)  +  Y*Y  +  (XALT-Z) * (XALT-Z) ) 
FD1= ( (SEP-X) » ( VTX-VXX) +Y* (VXY-VTY) ♦ (XALT-Z>  t (VTZ-VXZ> ) 
FDl-  FDl/<BPtLA) 

LN»C/<F0*1E9+FD1) 

FD2  »“1*<VTX»X  +  VTY»Y  VTZ»  (Z-RALT)  >  /  (APtLN) 
6ATE(I,4)-FD1  ♦  FD2 
CONTINUE 

FDMAX=MAX (GATE (1,4) , GATE (2, 4) , GATE (3,4) ,6ATE (4,4)) 
FDMIN=MIN (GATE (1,4) , GATE (2, 4> , GATE (3, 4) , GATE (4,4) ) 
DOPSR=ABS (FDMAX-FDMIN) 

RETURN 

END 


This  routine  prints  out  a  table  of  values  on  the  screen 

The  values  are  the  doppler  frequency  of  the  plane  of  interest. 

A  doppler  map  is  construe vted  by  hand  drawing  the  isodoppler 
contours. 

SUBROUT I NE  DOPPLT ( TABLE  > 

DIMENSION  TABLE (300, 250) , Iy( 15) 
do  30  1=1,11 

ly (I)=-150+25»i 
continue 

Print9,('  Idodoppler  Inf or eat i cm* ) 

Print*,  (»<X  axisX  Y  axis 

C  >’) 

Print  100, (IY(i),i«l,ll) 
do  70  i-1, 300,5 
li-i-51 

Print  105,ii, (TABLE(I,J>,J-1,11) 
continue 

forest (’  ‘,1116) 

forest ( 18, 11F6.0) 

RETURN 

END 

The  routine  for  interactively  changing  the  variables  of  interest 

SUBROUTINE  CHANGE  (DUMMY,  VJ,  TAZ,  TEL,  VX,  JtA2,  JtEL,  JPULSE,PRF,»rRF^ 
cFO, SEP , Z , XALT, P, GT,  6R, SNR, RCS, BW, LOSS, TEMP, Ri«L , RWAZ , XWEL, 
CXWA2,RALT) 

INTEGER  PRF,NPRF,D 


lor 


REAL  LOSS 
11  D«INT (DUMMY) 

22  GO  T0(600, 700, 800, 1600, 500, 100, 200, 300, 1400,900,1500,1300, 

C  2500,2400, 1000, 1200, 1 100, 2100, 1800,2000,400, 1700, 1900,2300, 

C  2200, 2600)  D 

100  PRINT*, » INPUT  TRANSMITTER  VELOCITY  (0-1200)' 

140  READ*,VX 

150  IF  (VX  .LT.  0)  THEN 

155  CALL  ERROR 

160  GO  TO  100 

170  ELSE  IF  (VX  .6T.  1200)  THEN 

175  CALL  ERROR 

177  60  TO  100 

180  END  IF 

190  GO  TO  2700 

200  PRINT*, » INPUT  TRANSMITTER  AZIMUTH  ANGLE  (0-360  DEGREES)' 

240  READ*,XAZ 

250  IF  (XAZ  .LT.  0)  THEN 

255  CALL  ERROR 

260  60  TO  200 

270  ELSE  IF  (XAZ  .GT.  360)  THEN 

275  CALL  ERROR 

277  GO  TO  200 

280  END  IF 

290  GO  TO  2700 

300  PRINT*, ’ INPUT  TRANSMITTER  ELEVATION  ANGLE  (-10  TO  +10  DEGREES)' 

340  READ1 , XEL 

350  IF  (XEL+IO  .LT.  0)  THEN 

355  CALL  ERROR 

360  GO  TO  300 

370  ELSE  IF  (XEL  .GT.  10)  THEN 

375  CALL  ERROR 

377  60  TO  300 

380  END  IF 

390  GO  TO  2700 

400  PRINT*, ' INPUT  TRANSMITTER  RECEIVER  SEPARATION  (10-200  KM)' 

440  READ*, SEP 

450  IF  (SEP  .LT.  10  )  THEN 

455  CALL  ERROR 

460  GO  TO  400 

470  ELSE  IF  (SEP  .GT.  200)  THEN 

475  CALL  ERROR 

477  GO  TO  400 

480  END  IF 

f'r')  GO  TO  2700 

0  PRINT*,  MNPUT  ALTITUDE  PLANE  OF  INTEREST  (0-30  KM)' 

540  READ*,Z 

550  IF  (Z  .LT.  0  )  THEN 

^55 - -CALL  ERWm - - - - — . . . . . . . . 

560  GO  TO  500 

570  ELSE  IF  (2  .GT.  30)  THEN 

575  CALL  ERROR 
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577  GO  TO  500 

580  END  IF 
590  GO  TO  2700 

600  PRINT*, 'INPUT  TARGET  VELOCITY  (0-2500)' 

640  READt,VT 

650  IF  (VT  .LT.  0)  THEN 

655  CALL  ERROR 

660  GO  TO  600 

670  ELSE  IF  (VT  .6T.  2500)  THEN 

675  CALL  ERROR 

677  GO  TO  600 

680  END  IF 

695  GO  TO  2700 

700  PRINT*,' INPUT  TARGET  AZIMUTH  ANGLE  (0-360  DEGREES)' 

740  READ* , TAZ 

750  IF  (TAZ  .LT.  0)  THEN 

755  CALL  ERROR 

760  GO  TO  700 

770  ELSE  IF  (TAZ  .GT.  360)  THEN 

775  CALL  ERROR 

777  GO  TO  700 

780  END  IF 

795  60  TO  2700 

800  PRINT*,' INPUT  TARGET  ELEVATION  ANGLE  (-10  TO  +10  DEGREES) 

840  READ* , TEL 

850  IF  (TEL+10  .LT.  0)  THEN 

855  CALL  ERROR 

860  GO  TO  800 

870  ELSE  IF  (TEL  .GT.  10)  THEN 

875  CALL  ERROR 

877  GO  TO  800 

880  END  IF 

890  GO  TO  2700 

900  PRINT*, » INPUT  TRANSMITTER  ALTITUDE  (0.01-30  KM)' 

940  READ* , XALT 

950  IF  (XALT  .LT.  0.01  )  THEN 

955  CALL  ERROR 

960  GO  TO  900 

970  ELSE  IF  (XALT  .GT.  30)  THEN 

975  CALL  ERROR 

977  60  TO  900 

980  END  IF 

990  GO  TO  2700 

1000  PRINT*,' INPUT  PU.SE  REPETITION  FREQUENCY  (1-50000)' 

1040  READ*,PRF 
1050  IF  <PRF  .LT.  1)  THEN 
1055  CALL  ERROR 

1060  60  TO  lOOO 

-  -1070  ^0000)  THEN - -  -  - 

1075  CALL  ERROR 

1077  60  TO  1000 

1080  END  IF 
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1090 
1100 
1140 
1150 
1155 
1160 
1170 
1175 
1177 
1180 
1190 
1200 
1240 
1250 
1255 
1260 
1270 
1275 
1277 
1280 
1290 
1300 
1340 
1350 
1355 
1360 
1370 
1375 
1377 
1380 
1390 
1400 
1440 
1450 
1455 
1460 
1470 
1475 
1477 
1480 
1490 
1500 
1540 
1550 
1555 
1560 
1570 
1575 
-a577  - 
1580 
1590 
1600 


GO  TO  2700 

PRINT#,  *  INPUT  NUMBER  OF  PRF  FREQUENCYS  <1-10)' 

READ»,NPRF 

IF  (NPRF  .LT.  1)  THEN 
CALL  ERROR 
GO  TO  1100 

ELSE  IF  (NPRF  .GT.  10)  THEN 
CALL  ERROR 
GO  TO  1100 

END  IF 

GO  TO  2700 

PRINT#, 'INPUT  PULSE  WIDTH  (0.1  TO  30  MICROSECONDS)' 

READ#,TPULSE 

IF  (TPULSE  .LT.  0.1)  THEN 
CALL  ERROR 
GO  TO  1200 

ELSE  IF  (TPULSE  .GT.  30)  THEN 
CALL  ERROR 
60  TO  1200 

END  IF 

60  TO  2700 

PRINT#, ’INPUT  CARRIER  FREQUENCY  (S  BAND  2-4  GHZ)’ 

READ#,F0 

IF  (FO  .LT.  2)  THEN 
CALL  ERROR 
GO  TO  1300 

ELSE  IF  (FO  .GT.  4)  THEN 
CALL  ERROR 
GO  TO  1300 

END  IF 

GO  TO  2700 

PRINT*,  ’  input  TRAI4SHITTER  POWER  (0.00050-2.0  MEGAWATTS)' 

READ#,P 

IF  (P  .LT,  0.00050  >  THEN 
CALL  ERROR 
GO  TO  1400 

ELSE  IF  (P  .GT.  2.0)  THEN 
CALL  ERROR 
GO  TO  1400 

END  IF 

60  TO  2700 

PRINT#, 'INPUT  TRANSMITTER  GAIN  (0-45  dB)' 

REA0#,6T 

IF  (GT  .LT.  0  )  THEN 
CALL  ERROR 
60  TO  1500 

ELSE  IF  (GT  .GT.  45)  THEN 
CALL  ERROR 

- 60  -JO  1500 - - - 

END  IF 

60  TO  2700 

PRINT#,  '  INPUT  MINIMUM  TAfcGET  CROSS  SECTION  (0.1-100  Sq  M)  ’ 
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1640  READ>,RC8 

1650  IF  (RC8  .LT.  0.1)  THEN 

1655  CALL  ERROR 

1660  GO  TO  1600 

1670  ELSE  IF  (RCS  .GT.  100)  THEN 

1675  CALL  ERROR 

1677  GO  TO  1600 

1680  END  IF 

1695  60  TO  2700 

1700  PRINT*, ’INPUT  RECEIVER  ANTENNA  GAIN  (0-45  dB) ' 

1740  READ*,6R 

1750  IF  (GR  .LT.  0)  THEN 

1755  CALL  ERROR 

1760  GO  TO  1700 

1770  ELSE  IF  (GR  .GT.  45)  THEN 

1775  CALL  ERROR 

1777  GO  TO  1700 

1780  END  IF 

1795  60  TO  2700 

1800  PRINT*, *  INPUT  RECEIVER  MINIMUM  DETECTION  SNR  (1-20  dB) ’ 

1840  READ*, SNR 

1850  IF  (SNR  .LT.  1)  THEN 

1855  CALL  ERROR 

1860  60  TO  1800 

1870  ELSE  IF  (SNR  .GT.  20)  THEN 

1875  CALL  ERROR 

1877  GO  TO  1800 

1880  END  IF 

1890  GO  TO  2700 

1900  PRINT*, ’ INPUT  RECEIVER  BANDWIDTH  (0.5-2.O  fCGAHERTZ) * 

1940  READ*,BW 

1950  IF  (BW  .LT.  0.5  )  THEN 

1955  CALL  ERROR 

1960  60  TO  1900 

1970  ELSE  IF  (BW  .GT.  2,0)  THEN 

1975  CALL  ERROR 

1977  GO  TO  1900 

1980  END  IF 

1990  GO  TO  2700 

2000  PRINT*,  •INPUT  ESTIMATE  OF  SYSTEM  LOSSES  (1-25  dB)' 

2040  READ*, LOSS 

2050  IF  (L(^  .LT.  1)  THEN 

2055  CALL  ERROR 

2060  GO  TO  2000 

2070  ELSE  IF  (LOSS  .GT.  25)  T)CN 

2075  CALL  EKFiOR 

2077  GO  TO  2000 

2080  END  IF 

2090  -GO  TO  2700 - -  -  - - - - -  “  '  - - - 

2100  PRINT*, 'INPUT  RECEIVER  FRONTEND  TEMPERATURE  (290-500 
CDEGREES  K)  ’ 

2140  READ!, TEMP 
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2150 
2155 
2160 
2170 
2175 
2177 
2130 
2190 
2200 
2240 
2250 
2255 
2260 
2270 
2275 
2277 
2280 
2295 
2300 
2340 
2350 
2355 
2360 
2370 
2375 
2377 
2380 
2395 
2400 
2440 
2450 
2455 
2460 
2470 
2475 
2477 
2480 
2495 
2500 
2540 
2550 
2555 
2560 
2570 
2575 
2577 
2580 
2593 
.  2600 
2605 
2610 
2615 


IF  (TEMP  .LT.  290  TKEN 
CALL  ERROR 
GO  TO  2100 

ELSE  IF  (TEMP  .6T.  500)  THE>I 
CALL  ERROR 
GO  TO  2100 

END  IF 

60  TO  2700 

PRINTt,* INPUT  RECEIVER  ANTENNA  BEAM  HEIGHT  (1-8  DEG)’ 
READ*tRWEL 

IF  (RWEL  .LT.  1)  THEN 
CALL  ERROR 
GO  TO  2200 

ELSE  IF  (RWEL  .GT.  8)  THEN 
CALL  ERROR 
GC  TO  2200 
END  IF 
60  TO  2700 

PRINT*,  MNPUT  RECEIVER  ANTENNA  BEAM  WIDTX  (1-5  DEG)’ 
REAO«,RWAZ 

IF  (RWAZ  .LT.  1)  THEN 
CALL  ERROR 
GO  TO  2300 

ELSE  IF  (RWAZ  .GT.  5)  THEN 
CALL  ERROR 
60  TO  2300 
END  IF 
60  TO  2700 

PRINT*, *  INPUT  TRANHITTER  ANTENNA  BEAM  HEIGHT  (1-8  DEG) 
READ* , XWEL 

If  (XWEt.  .LT,  1)  THEN 
CALL  ERROR 
GO  TO  2400 

ELSE  IF  (XWEL  .GT,  8)  THEN 
CALL  ERROR 
60  TO  2400 
END  IF 
60  TO  2700 

PRINT*,  *  INPUT  TRANHITTER  ANTEMM  BEM1  WIDTH  <1-5  DE&)’ 
REAOt.XWAZ 

IF  (XWAZ  .LT.  1)  THEN 
CALL  ERROR 
GO  TO  2500 

ELSE  IF  (XWAZ  .GT.  5)  THEN 
CALL  EW^OR 
GO  TO  2500 
END  IF 
GO  TO  2700 

PRINT*, ’INPUT  RECEIVER  JiLilTUDE  10. Oi  30  KlLOIffiTERS) 

READ*,  RALT 

IF  (RALT. LT, 0,01)  THEN 
CALL  ERROR 


2620  GO  TO  2600 

2625  ELSE  XF  (RALT.GT.30>  THEN 

2630  CALL  ERROR 

2635  GO  TO  2600 

2640  END  IF 

2700  RETURN 

2210  END 

c 

c  the  routine  to  calculate  the  range  constant  4or  the  power 
c  limiting  equation  region 
c 

1  SUBROUTINE  RANGE <P,6T,6R, SNR, RCS,F0, LOSS. 6W, TEMP, RAN6E4) 

10  REAL  C2, LOSS, TEMP, GT,GR,P,RCS,F0,BW,RAN6E4 

15  C2»(6R+6T~L0SS-SNR)/10 

20  C2=10»*C2 

25  C3=3. 2865E6 

30  RANGE4=P»RCS»C2tC3/ (BWtTEMPtFOtFO) 

35  RETURN 

40  END 

c 

c  The  :'outine  to  cacluation  the  doppler  frequencies  for  the  plane 
c  of  in«'©rest 

c 

1  SUBROUTINE  i ORMAP <VX, XEL, XAZ , VT, TEL, TA2 , F0,C, SEP. XALT, 

CZ, TABLE, RALT) 

5  DIMENSION  TABLE '300,2501 

10  REAL  LA.LN 

600  VXZ  =  VX «S I ND(XEL 1/3600 

6C5  VXY  =  VXtCOSD<XELltSIND(XAZ 1/3600 

610  VXX  =  VX«C0SD(XELltC0SD<XA2 1/3600 

615  VTZ  «  VT »S I ND( TEL 1/3600 

620  VTY  •=  VTtCOSD (TEL  1*SIN0<TA2 1/3600 

625  VTX  =  VT »COSD{ TEL )»C0SD(TA2 1/3600 

630  LA«C/(F0*IE91 

o35  DO  710  1^1,300,5 

640  X-I~51 

645  00  70S  J»1,U 

650  V— t50^jt25 

655  fiP  -  SQRTfXiX  ♦  Y»Y  ♦  (Z-RALT) • IZ-RALT) 1 

660  BP  •  SORT  ((SEP~X>9(SEP<-X)  ♦  VtV  ♦  tXALT~Z)«(XALT-Z) ) 

674  FDi«((SEP-Xl<(VTX''VXX)^Y$<VXY^TY>a<XALT-Z>t(VXZ~VXZ)> 

675  FDl«FDl/(B«»tLAl 

677  LN*C/ (FO* lEV^FDl > 

680  FD2  =-!•  tVTXtX  ♦  VTVtY  ♦  VTZt (Z-RALTl )  /  <AP*LN) 

700  TABLE(I,J1-FD1  ♦  FD2 

705  CONTINLC 

710  CONTINUE 
735  RETURN 

740  -END  - -  - - - - - -  — - 

c 

c  The  routine  that  calculates  the  minimum  and  max i mum  values 
c  for  the  isorange  ellipses 


10*? 


c 

1  SUBROUTINE  MlNilAX (TPUL8£»C,8EP,NPRF,PRF,  180NIN,  ISOHAX,  I, 

C  XALT,RALT) 

S  REAL  ISOMIN, X80NAX 

7  i nteger  PRF 

e  R3  «  SQRKSEPtSEP  ♦  (XALT-RALT)  tt2) 

10  IS0MIN*2*TPULSE*lE-6tC^R3 

15  ISOMAX- (NPRFtC/PRF) ♦RS- (TPULSE* 1E-A$C) 

20  I-INT( (IS0MAX-IS0MIN)/10)+2 

25  RETURN 

30  END 
C 

C  This  routine  plots  the  ov«ls  of  Cessini 

c 

1  SUBROUTINE  POWER(2,XALT,SEP,R4,BASE,RALT> 

5  DIMENSION  UPPER (302), DOWNER (302), BASE (302) 

6  UPPER(1)»301 

7  DOWNER ( 1 ) -301 

11  C4»SEP»(-2) 

DO  135  1—50,250 
X-I 

31  C1«(Z-RALT)I»2 

72  C2=(Z-XALT)»*2 

73  C3=C2+Cl+SEP4t2 

741  C5-C4$C1 

742  C6«CltSEP$12+CltC2-R4 

122  B4-X»»2^C3/2-X$SEP 

123  CC-Xtt4«C4tX«t3<*^C3»X»42^CS«X^C6 
B5-B4t»2-CC 

IF(CC.GT.O)  THEN 
UPPER (I ♦52) «0 
GO  TO  130 
end  if 

124  UPPER (X^52) -SORT (-B4 ♦SORT (B5)) 

130  Dawr«lER(I^52)— UPPER  (I^52) 

135  CONTINUE 

232  CALL  CPLOT( BASE. UPPER) 

233  CALL  CPLOT (BASE, DOWNER) 

235  RETURN 

240  END 
c 

c  This  routine  plots  the  tsorenge  ellipsiotfs  snd  keeps  treck 
c  of  n^ethcr  or  not  the  lest  required  ellipsiod  hes  been  tfreen 

c 

1  SUBROUTINE  ISORdSOMIN.  ISOHAX,  I.SEP,  Z,XALT,6AS£,RALT) 

5  DIMENSION  UPPER (302) .DOWNER (302) .BASE (302) 

6  REAL  ISOMIN, ISOMAX 

7  UPPER(1)-301 

..B  .  DOWNER  (1)  *301  _ _ _ _ _  . . - 

10  SEP2-SeP/2 

GO  TO { 100, 200, 300, 300, 300, 300, 3(0, 300, 300, 300, 300, 300, 
C300,  300,  300, 300, 300,  300,  300, 300, 300,300,300, 300,  300, 
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C300,300,300,300»300,300,300)  I 
SOO  CURVE-I80M1N 
60  TO  400 
200  CURVE^ISOMAX 
60  TO  400 

300  CURVE«10#<I-2*INT(ISOHIN/J0)) 

400  I«I-1 

20  A*CURVE/2 

23  A2«A«A 
B2=A2~SEP2tSEP2 
C12-B2/A2 

24  sTHETA= ( XALT-RALT) /SEP 

25  CTHETA=COS (ASIN (THETA) > 

26  ClC«Z*CTHETA-RALTtCT)CTA 

27  Cll=Z*STHETA-RALT«STHEtA>fSEP/2 

45  DO  65  J=-50,250 

46  X«=J 

50  Y=  (B2-  (XtSTHETA-^ClO)  •»2-Cl2»  (XtCTHETA-Cl  1 )  tt2> 

i-f  (y.lt.O)  then 
upper ( j *52) “0 
go  to  55 
end  if 

upper  ( j+52)  ^^sgrt  (y ) 

55  downer  ( j+52)  =~upper  ( j-*'52) 

65  CONTINUE 

66  CALL  Lined +2) 

67  CALL  CPLOT< BASE, UPPER) 

68  CALL  CPLOT(BASE,DOWf«:R) 

70  RETURN 

75  END 

c 

c  This  routine  clears  the  40XX  screen  when  in  the  40XX  mode, 
c 

1  SUBRCMJT I  HE  WIPE 

5  CALL  INITT<1200) 

10  CALL  BINITT 

15  CALL  FINITT (0,700) 

20  RETURN 

25  END 
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